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ABSTRACT 

We present the first statistical study of X-ray cavities in distant clusters of galaxies 
{z > 0.3). With the aim of providing further insight into how AGN feedback operates 
at higher redshift, we have analysed the Chandra X-ray observations of the Massive 
Cluster Survey (MACS) and searched for surface-brightness depressions associated 
with the Brightest Cluster Galaxy (BCG). The MACS sample consists of the most 
X-ray luminous clusters within 0.3 ^ ^ ^ 0.7 (median Lx,rass = 7 x 10^^ erg s~^ ), 
and out of 76 clusters, we find 13 with "clear" cavities and 7 with "potential" cavities 
(detection rate ~ 25 per cent). Most of the clusters in which we find cavities have a 
short central cooling time below 3 — 5 Gyrs, consistent with the idea that cavities sit 
predominantly in cool core clusters. We also find no evidence for evolution in any of the 
cavity properties with redshift, up to z ~ 0.6. The cavities of powerful outbursts are 
not larger (or smaller) at higher redshift, and are not able to rise to further (or lesser) 
distances from the nucleus. The energetics of these outbursts also remain the same. 
This suggests that extreme "radio mode" feedback (Lmech > 10^"^ erg s~^ ) starts to 
operate as early as 7 — 8 Gyrs after the Big Bang and shows no sign of evolution since 
then. In other words, AGNs lying at the centre of clusters are able to operate at early 
times with extreme mechanical powers, and have been operating in such a way for at 
least the past 5 Gyrs. 

Key words: X-rays: galaxies: clusters - cooling fiows - galaxies: jets 



1 INTRODUCTION 

Active galactic nucleus (AGN) feedback plays a major 
role in quenching star formation, enriching the surrounding 
medium with metals, and fuelling the supermassive black 
hole (SMBH) of the host galaxy. Some of the most extreme 
examples of AGN feedback are seen in clusters of galaxies, 
where the central SMBH inflates large cavities filled with 
radio emitting particles through jets. These cavities appear 
as depressions in the X-ray image and provide a direct mea- 
surement of the energy being injected into the surrounding 
medium by the central SMBH (see reviews on the topic by 
Peterson & Fabian 2006 and McNamara & Nulsen 2007) . 
Detailed studies of individual nearby clusters such 
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as the Perseus Cluster (see iFabian et al. 1 l2006l . and 
references therein) have shown that the AGN ly- 
ing at the centre can energetically offset cooling of 
the intracluster medium (ICM) not only by inflating 
these large cavities, but also by inducing weak shocks 
and propagating e n ergy through sound /pressure waves 



(iBirzan e t all l2004l: iRaffertv et al.1 l2006l: IPunn FabianI 



2006, 200 81: iMcNamara NulsenI l2007l: iFabian et al.ll2003l . 
.2006 : .Forman et al.ll2005l : ISanders Fabianll2007l l . Statisti- 
cal studies of X-ray cavities have also provide d a wealth of 
information on how AGN feedback o perates (Birz an et al.l 



2004. 200 8: 'Punn et al.1 l2005l: jPunn Fabian 20oi 



Nulsen et al. 2007; Punn et al 



2010y Cavagnolo et al.l 



200i 



201 



Pong et al.. .20ia : .O'Sullivan et al.ll201lh . Essentially, these 
studies find that the power output of cavities is substantial, 
in most cases sufficient to prevent cooling of the ICM, and 
that cavities are common, especially in systems which re- 
quire some form of heating to prevent the gas from cooling. 
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However, all the statistical compilations of cavities have un- 
til now focussed on the nearby Universe (only 4 objects are 
at 2; > 0.3). The cavity properties at the higher redshift end 
therefore remain unexplored. It is also not clear if and how 
AGN heating evolves acros s time . According to cosmo logical 
simulations (|Croton et al.1 12OO6I : ISiiacki Springell I2OO6I ). 
present-day black holes are in a "radio mode" phase, which 
involves a black hole accreting at sub-Eddington rates and 
driving powerful outflows, whereas at earlier times ^ 1), 
"quasar mode" feedback dominates and consists of a merger- 
driven phase where the black hole grows rapidly. When and 
how this transition occurs remains unclear. 

We can start investigating how AGN feedback evolves 
with redshift by using targeted samples oi z > 0.3 clus- 
ters that have extensive follow-up X-ray observations, and 
in this case, the Massive Cluster Survey (MACS) provides 
an ideal example. MACS was launched in 1999 and compiled 
the first large sample of very X-ray luminous , and therefore 
very massive, clus ters in the distant Universe (|Ebeling et al.l 
I2OOIL I2OOI I2OIQI ). The MACS sample was based on the 
ROSAT All-Sky Survey (RASS; ffriimper 1983, and refer- 
ences therein) Bright Source Catalogue (BSC; Voges et al. 
1999). It consists of clusters above a redshift of 2; > 0.3 
and covers a solid angle over 22000 deg^ in the extra- 
galactic sky. In Fig. [U we show the Lx — z distribution 
for various samples, includ ing the Brightest Cluster Survey 
(jEbeling et al.lll998l. I2OOOI). the Einstein M edium Sensitiv- 



eBCS 



MACS 



ity Survey (EMSS; ICioia k Luppinol 



19941). the Wide An- 



gle Rosat Pointed Survey f WARPS: iPerlman et aDl2002l ). 
the 400 deg^ cluster samples (Burenin et al. 2007), and the 
MACS sample. This figure shows that MACS selects on av- 
erage clusters which are 10 — 20 times more luminous than 
the other samples in the same redshift range. The MACS 
sample now contains 124 spectroscopically confirmed clus- 
ters within 0.3 ^ ^ ^ 0.7, for which more than two thirds 
are new discoveries. 

Since MACS contains only very X-ray luminous clus- 
ters (median Lx,rass = 7 x 10^^ erg s~^ ), those harbour- 
ing a cool core should require extreme feedback from their 
central AGN to prevent the ICM from cooling (Lmech ^ 
-|^q44_45 gj.gg-i^_ therefore expect to see some of the 
most extreme examples of AGN outbursts in this sample, 
which will allow us to extend our knowledge of X-ray cavi- 
ties not only in the distant Universe {z > 0.3), but also at 
the higher end of mechanical power arising from the AGN. 
To date, more than two thirds of all MACS clusters have 
been targeted in Chandra follow-up observations, and the 
unique and exquisite spatial resolution of Chandra is essen- 
tial to identify cavities at such redshifts. The MACS sample 
therefore provides a unique opportunity to study the cav- 
ity properties in clusters at z > 0.3, and determine if AGN 
feedback evolves with time. 

We first present the data in Section 2, and then discuss 
how we proceeded in identifying systems with cavities in Sec- 
tion 3. Section 4 and Section 5 outline the techniques used 
to calculate the energetics of the cavities, as well as their as- 
sociated time-scales. In Section 6, we derive the properties 
of the cores and plot scaling relations in Section 7. Finally, 
we present the results in Section 8 and discuss them in Sec- 
tion 9. We adopt Hq = 70 km s~^ Mpc"^ with Qm = 0.30 
and = 0.7 throughout this paper. All errors are 2cr, and 



^ 10.0- 





0.3 0.4 

redshift z 



Figure 1. The X-ray luminosity versus redshift distribution of 
different cluster samples. There are a total of 124 MACS clusters, 
all with z > 0.3. The MACS clusters outlined with black circles 
are those in which we identified cavities. 



the abundance ratios of lAnders k Grevessd (| 19891 ) were used 
throughout this paper. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 X-ray observations 

The X-ray data were obtained through the Chandra data 
archive. We found observations with the Advanced CCD 
Imaging Spectrometer (ACIS) for 83 of the MACS clus- 
ters, with exposure times varying between 10 ks and 100 ks. 
Some objects were also observed multiple times. For each 
cluster, we selected the deepest observational data set 
available. These data were then processed, cleaned and 
calibrated using the latest version of the CIAO software 
(ciAOv4.3, CALDB4.4.1), and starting from the level 1 event 
file. We applied both CTI (charge time interval) and time- 
dependent gain corrections, as well as removed flares using 
the LC_CLEAN script, with a 3a threshold. When a cluster 
was observed multiple times with the target centred on the 
same detector and in the same observing mode (FAINT or 
VFAINT), we combined the different observations only if 
this would improve the image quality significantly. We then 
exposure-corrected the images, using an exposure map gen- 
erated with a monoenergetic distribution of source photons 
at 1.5 keV( which is almost the peak energy of the clusters). 
These images were used to identify systems with cavities (see 
Section 3). Out of the 83 clusters for which we have X-ray 
observations, 2 are clear mergers with multiple components 
merging together and 5 have a bright point source with no 
clear extended emission associated with it. We therefore only 
consider the remaining 76 in the analysis of X-ray cavities. 

For the clusters in which we identified cavities (see Ta- 
ble [p, the spectr a were analysed using XSPEC (vl 2.6.0 d,f, 
e.g. Arnaudiri996l V For each object, the Galac tic absorption 
was kept frozen at the iKalberla et al. I EM) value. Vary- 
ing this parameter did not improve significantly the fit in 
any of our clusters. For targets observed multiple times, we 
fitted simultaneously the spectra. Finally, although some of 
the objects in our sample with cavities have a bright X-ray 
point source at their centres, we find no significant pileup 
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Table 1. Chandra observations of the MACS clusters with cavities - (1) Name; (2) Alternate name; (3) Redshift z; (4) Observation 
identification number; (5) Exposu re time; (6) Chip nam e w here the targe t was centre d on; (7) Refe rences where the redshift was taken 
from: (i) Ebeling et al. (2007); fii) lEbehng et al.l (l201Cj): fiii ^ lStocke et al.l (Il99lh : (W) IWright et al.l (1983); (v) Mann &; Ebeling ( 2m|); 
(vi) Kleinmann et al.l ((19881 ): (vii) iDressle^^Cunnl (j 19921 ): (viii) Ebeling et al. (in preparation); (ix) Sloan Digital Sky Survey; (x) 
[Smail et al.. (,2007 ). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Cluster 


Alternate 


z 


Obs. ID 


Exposure 


Chip 


Ref. 


name 


name 






ks 






MACS J0159.8-0849 




0.404 


3265 


14.5 


ACIS-I3 


(ii) 








6106 


28.9 


ACIS-I3 










9376 


18.0 


ACIS-I3 




MACS J0242. 5-2132 




0.314 


3266 


8.6 


ACIS-I3 


(iv) 


MACS J0429. 6-0253 




0.397 


3271 


21.4 


ACIS-I3 


(ii) 


IVLALyo JUo4/.U-oyU4 




U.oiy 


Q07Q 


ly.z 






MACS J0913. 7+4056 


IRAS 09104+4109 


0.442 


10445 


70.4 


ACIS-I3 


(vi) 


MACS J0947.2+7623 


RBS 0797 


0.354 


7902 


38.8 


ACIS-S3 


(ii) 


MACS J1411. 3+5212 


3C295 


0.460 


2254 


76.4 


ACIS-I3 


(vii) 


MACS J1423. 8+2404 




0.5449 


4195 


106.9 


ACIS-S3 


(i) 


MACS J1532. 8+3021 


RX J1532. 9+3021 


0.3613 


1649 


9.5 


ACIS-S3 


(ii) 


MACS J1720. 2+3536 


Z8201 


0.3913 


3280 


17.6 


ACIS-I3 


(ii) 








6107 


27.3 


ACIS-I3 










7718 


6.8 


ACIS-I3 




MACS J1931. 8-2634 




0.352 


9382 


95.0 


ACIS-I3 


(ii) 


MACS J2046. 0-3430 




0.423 


9377 


35.9 


ACIS-I3 


(viii) 


MACS J2140. 2-2339 


MS 2137.3-2353 


0.313 


4974 


7.3 


ACIS-S3 


(iii) 








5250 


34.7 


ACIS-S3 










928 


34.8 


ACIS-S3 




MACS JOlll. 5+0855 




0.485 


3256 


15.2 


ACIS-I3 


(v) 


MACS J0257.1-2325 




0.5039 


1654 


18.0 


ACIS-I3 


(i) 








3581 


16.1 


ACIS-I3 




MACS J1359. 1-1929 




0.447 


9378 


46.1 


ACIS-I3 


(v) 


MACS J1359.8+6231 


MS 1358.4+6245 


0.330 


516 


48.5 


ACIS-S3 


(iii) 


MACS J1447.4+0827 


RBS 1429 


0.3755 


10481 


11.5 


ACIS-S3 


(ix) 


MACS J2135. 2-0102 


IRXS J213515.7-010208 


0.325 


11710 


24.7 


ACIS-I3 


(x) 


MACS J2245. 0+2637 




0.301 


3287 


11.5 


ACIS-I3 


(ii) 



for any of these point sources. Pileup occurs when two or 
more photons are detected as o ne event (see for more details 
iDavisI 12001: 'Russell et alllioTo ). and the amount of pileup 
can be estimated by comparing the fraction of good (grades 
0,2,3,4,6) to bad grades (grades 1,5,7) for each point source. 
Typically, pileup becomes problematic when the fraction of 
bad grades exceeds 10 per cent of the good grades, but we 
find a ratio less then 6 per cent for all of our central point 
sources. 



2.2 Radio observations 

For each cluster with cavities, we searched for radio emission 
associated with the central AGN by looking through differ- 
ent radio surveys available to the public. We first looked 
for radio emission associated with the AGN at 1.4 GHz 
by using the 1.4 GHz VLA Faint Images of the Radio 
Sky at Twenty- Cent imeters surv ey (FIRST, average reso- 
lution of 5"; Becker et al.lll994l ). or the 1.4 GHz NRAO 
VL A Sky Survey catalo gue (NVSS, average resolution of 
45^^ ICondon et ahlllQQSl ). We then looked for radio emission 
associated with the AGN at lower frequency, either with 



the 326 MHz Westerb ork Northern Sky Sur ve.^B (WENSS, 
resolution of ^ 50^^ iRengelink eraP Il997l \ the 74 MHz 
VLA Low- frequency Sky Survey (VLSS, resolution of ^ 80''; 
ICohen et al. 2007), the 150 MHz TIER GMRT Sky Survey 
(TGS^l, resolution of - 20'') or the 843 MHz Sydney Uni- 
versity Molonglo Sky Survey (SUMS S, resolution of ^ 40"; 
iBock et al.lll999l : iMauch et al.|[2003h . Finally, we searched 
for high-frequency radio emission associated with the cen- 
tral AGN at 5 GHz with the Parkes-MIT-NRAO (PMN, 
resolution of ^ 5') radio survey (Grif fith Wrightlll993l l or 
the EINSTEIN Observatory Extended Medium- Sensitivity 
Survey (EMSS, wher e the radio resolut i on of the survey lie s 
between - 4 - 14"; iGioia et al.lll990l : IStocke et aP 1 199 il l, 
and then at 28.5 GHz with the Berkeley-Illinois-Maryland 
Association (BI MA, sensitive to angular scales around 1.5'; 
ICoble et al. I [2007. ). If no point source was seen within the 
central regions (r ^ 100 kpc), we used the 2crrms value within 
the beam area as an upper limit for the flux. The values we 
obtain are shown in Table [2] Some sources also have flux 
densities available from the literature, and we include these 
in Table 



^ http:/ /www.astron.nl/wow/testcode.php?survey=l 
^ http://tgss.ncra.tifr.res.in 
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Errors are derived as the quadratic sum of the rms 
noise level in the map and the systematic uncertainty as- 
sociated with the value. For some of the surveys used, the 
error quoted in the catalogue did not account for systematic 
uncertainties. This type of error va ries with frequency , but 
is on the order of 5 per cent fsee ICarilli et al.lll99ll ). For 
simplicity, we therefore choose to compute the total uncer- 
tainty assuming a 5 per cent systematic error and a 2crrms 
noise level. MACS J1931. 8-2634 has a complicated morphol- 
ogy at radio wavelengths, with a Narrow Angle Tail (NAT) 
source located Ab" to the south that is contaminating the 
flux measurements of the cent ral g!:alaxy. The 1.4 GHz flux 
quoted in Table [2] is taken from lEhlert et al ] (|201lh and only 
includes the contribution of the central galaxy. However, for 
the 5 GHz and 150 MHz data, we had to remove the contri- 
bution of the NAT source. In this case, we used the peak flux 
value of the galaxy and subtracted the peak contribution of 
the NAT source. 

Finally, by using the flux densities quoted in Table (2] 
we calculated a rough estimate of the total radio luminosity 
(i^radio) with Eq. [T] and integrating between — 10 MHz 
and z/2 = 10000 MHz. 

Lradio = 47rL>£ / \s.)dv (l) 

Here, Dl is the luminosity distance to the source and Su 
is the flux density {S^ oc where a is the spectral in- 

dex). For the central galaxies of MACS J0242. 5-2132 and 
MACS J1411.3+5212 (3C 295), we also use the extensive 
follow-up radio observations available in the NASA/IPAC 
Extragalactic Database (NED) to determine Lradio- Only 
MACS J141 1.3+5212 has a complete coverage from 10 
MHz to 10 GHz. For the remaining objects, we extrapolate 
the values at 10 MHz and 10 GHz based on the assump- 
tion that the flux density scales as Su oc and for each 
extremity (10 MHz and 10 GHz), we use the local spectral 
index as determined from the two nearest flux density data 
points. Using a simple trapezoid rule, we then integrate over 
the range 10 MHz to 10 GHz using the extrapolated values 
at 10 MHz and 10 GHz, as well as the other data points 
available for each source. Our results are shown in Tabled 



2.3 Optical observations 

We also searched for optical observations of each clus- 
ter harbouring a cavity. These images were used to verify 
that the cavities identified were associated with the central 
dominant galaxy. For each cluster, we searched for Hub- 
ble Space Telescope (HST) archival images (many of which 
were acquired through snapshot campaigns, see Ebeling 
I2OO9I I, an d if not availa b le, w e used the R passband im- 
ages from ' Ebeling et all (|2007l . 120101 ) taken with the Uni- 
versity of Hawaii (UH) 2.2-m telescope (MACS J2046.0-3430 
and MACS JOll 1.5+0855) or the red images from t he Sec- 
ond P alomar Observatory Sky Survey [POSSII] Rei d et al] 
Il99lh for MACS J1447.4+0827. In Fig.H we show these op- 
tical images along with the 0.5 — 7.0 keV Chandra X-ray 
images. 



3 IDENTIFYING SYSTEMS WITH CAVITIES 

For each cluster, we visually searched for circular or ellip- 
soidal surface brightness depressions in the X-ray images 
associated with the central regions (r < 100 kpc). Beyond 
100 kpc, the count rate in most observations is not suffi- 
cient to identify cavities. As a first indicator, we computed 
unsharp-masked images for each source. This technique en- 
hances deviations in the original image and consists of sub- 
tracting a strongly smoothed image from a lightly smoothed 
image, and is similar to subtracting an elliptical model of 
the cluster emission. Both the unsharp-masked and ellipse- 
subtracted images are shown in Fig. lAll However, although 
these figures can be very useful, they can also be mislead- 
ing and cause false cavities to appear. Therefore, we only 
considered a cavity as being truly there if we could also 
see a hint of a depression in the original image. Once the 
systems with cavities were identified, we proceeded in clas- 
sifying each cavity (see Column 2 of Table (3]). First, if the 
cavity had a clear contrast in the raw image, as well as in the 
unsharp-masked and ellipse-subtracted image, we classified 
it as being "clear". Out of these "clear" cavities, some had 
surrounding bright rims in the raw image which we iden- 
tified with the annotation 1 in Table [3] Those without the 
bright rims were given the annotation 2. If there was just a 
hint of a depression in the raw image, but the cavity could 
clearly be seen in the unsharp-masked or ellipse-subtracted 
image, then we classified it as being a "potential" one. 

In total, we find that 13 clusters have at least one 
"clear" cavity and 7 have "potential" cavities. Therefore, 
out of the 76 MACS clusters with Chandra data, we find 
that 20 have X-ray cavities. If we calculate the fractional 
difference between the counts within the cavity and the sur- 
rounding region, then we find that the decrements are on 
the order of 20 — 30 per cent for the "clear" cavities and 
between 10 — 20 per cent for the "potential" cavities. These 
fractions are roughly consistent with the decrements seen in 
nearby clusters of galaxies. In the calculation, we use the 
two immediate regions opposite the cavity, at the same ra- 
dius and same size as the cavity, as the surrounding region. 
Note also that many of the MACS clusters have only been 
observed for 10 — 20 ks and some have as little as 500 — 1000 
counts within 200 kpc. We could therefore be missing many 
cavities and our detection rate should only be considered as 
a lower limit to the number of MACS clusters with cavities. 

Our sample of objects with cavities includes MACS 
1411.3+5212 (3C 295). This is a particularly interesting 
cluster, which has a very powerful radio source at its cen- 
tre and shows strong evidence for Inverse Compton scat- 
tered Cosmic Microwave Background (ICCMB) photons 
along the central jet axis. Out of the remaining clusters, 
MACS J2135. 2-0102 has a central galaxy harbouring a dou- 
ble or triple nucleus (see optical image in Fig. [2]). Initially, 
we identified two cavities in this system, one to the north- 
west and another to the south-east. However, the BCG is 
located almost within the south-eastern cavity. Although 
this could be due to projection effects, we choose to dis- 
card the southern cavity and just consider the northern one 
as a "potential" cavity. MACS J2135. 2-0102 is more weh 
known for its lensing a rcs in the shape of a "Cosmic Eye" 
and "Cosmic Eyelash " (|Smail et al.ll2007l : ISiana et al.ll2009l : 
ISwinbank eraL(l2010h . 
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Table 2. Radio properties of the MACS clusters with cavities - (1) Name; (2) 74 MHz radio flux density from VLSS; (3) 326 MHz 
radio flux density from WENSS; (4) 1400 MHz radio flux density; (5) 5 GHz radio flux density; (6) 28.5 GHz radio flux density from 
Berkeley-IUinois-Maryland Association (BIMA; Coble et al. 2007); (7) Other radio flux density; (8) Integrated radio luminosity from 10 
MHz to 10 GHz, see Section 2.2 for details; (9) References: (i) 1.4 GHz First; (ii) 1.4 GHz NVSS; (iii) Parkes MIT-NRAO 4.85 GHz 
survey fCriflith & Wright 1993); (iv) 5 GHz private communication (M. T. Hogan); (v) 5 GHz from Hines & Wills (199^); (vi) 5 GHz 
from Cavagnolo et al. ( 20111): fvii) 5 GHz fromlCregorv CondonI (Il991h : (viii) EMSS 5 GHz (|Gioia et al.lll990l : IStocke etaD [l991): (ix) 
843 MHz SUMSS ( Bock et a .lll999l:[Mauch et al.ll2003l ): fx) 15 GHz Arcminute Microkelvin Imager (AMI), private communication (K. 
Grainge); (xi) 1.4 GHz from lshlert et al.l (l201ll ): (^ii) 150 MHz TIFR GMRT Sky Survey "For the central galaxies of MACS J0242.5- 
2132 and MACS J1411. 3+5212 (3C 295), we used the extensive follow-up radio observations available in the NASA/IPAC Extragalactic 
Database (NED) to determine Lradio- ^MACS J1931. 8-2634 has a complicated morphology at radio wavelengths, see Section 2.2 for 
details on how we obtained the flux densities. 















\' ) 






Cluster 


5'74MHz 


'S'326MHz 






5'28.5GHz 


'S'other 


-^radio 


Ref. 


name 


mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


(10^2 erg s-i) 




MACS J0159.8-0849 


< 112 




31.4+1.6 


58+11 






3.56 


(i,iii) 


MACS 10242.5-2132'^ 


890 ± 145 




1255 + 73 


795 + 43 






27.1 


(ii,iii) 


MACS J0429.6-0253 


< 214 




138.8 + 8.1 








7.2 


(ii) 


MACS J0547.0-3904 






31.4+1.9 


15.4 + 0.8 




19.6 + 1.3[@843MHz] 


0.59 


(ii,iv,ix) 


MACS J0913. 7+4056 


< 272 


54.0 + 7.1 


8.3 + 0.5 


1.6 + 0.1 


0.69 + 0.12 


0.80 + 0.04[@15GHz] 


1.33 


(i,v,x) 


MACS J0947.2+7623 


< 329 


91.0 + 6.1 


21.7+1.3 


4.0 + 0.2 






1.02 


(ii,vi) 


MACS 11411.3+5212'^ 


120270 ± 6022 


61647 + 3082 


22171 + 1109 


7401 + 808 






1025.3 


(i,vii) 


MACS J1423.8+2404 


< 232 




5.2 + 0.4 




1.49 + 0.12 




3.58 


(i) 


MACS J1532.8+3021 


< 222 


71.0 + 8.2 


17.1 + 0.9 


8.8 + 0.5 


3.25 + 0.18 




0.94 


(i,iv) 


MACS J1720.2+3536 


< 266 


103.0 + 7.4 


16.8 + 1.0 








1.23 


(i) 


MACS 11931.8-2634^ 






70 + 4 


6.0+1.3 




2799 + 161 [@150MHz] 


80.4 


(iv,xi,xii) 


MACS J2046.0-3430 






8.1 + 0.6 






13 + 1.3[@843MHz] 


2.49 


(ii,ix) 


MACS J2140.2-2339 


< 116 




3.8 + 0.5 


1.0 + 0.1 






0.42 


(ii,viii) 


MACS JOlll. 5+0855 


< 124 




< 0.25 








< 3.20 


(i) 


MACS J0257.1-2325 


< 100 




< 0.89 








< 1.50 


(ii) 


MACS J1359. 1-1929 


< 266 




9.5 + 0.8 








2.34 


(ii) 


MACS J1359.8+6231 




< 7 


2.8 + 0.6 


3.8 + 0.2 


1.67 + 0.08 




0.16 


(i,viii) 


MACS J1447.4+0827 


< 203 




39.0 + 2.0 








2.00 


(i) 


MACS J2135. 2-0102 


< 192 




< 0.40 








< 1.90 


(i) 


MACS J2245.0+2637 


< 208 




5.8 + 0.6 








0.72 


(ii) 



Table 3. Cavity energetics. The annotations in Column 2 indicate the quality of the cavities detected. Those indicated with a ^ are the 
cavities with clear surface brightness depressions. The second annotation -"^'-^ indicates whether a cavity is surrounding by bright rims (^) 
or not (^). Column 3 shows the position angle of the major axis along the cavity, and is measured counter-clockwise from the north. 



Cluster Lobe PA Ri i?w R £^bubble = PV tc^ tbuoy ^refill ^'cs ^buoy -Prefill 

(degrees) kpc kpc kpc 10^^ erg 10^ yrs 10^ yrs 10^ yrs IC^^ erg s"-*^ lO'^'^ erg s"-*- IC^^ erg s"-*^ 



MACS J0159.8-0849 




29 


16.4 


8.8 


25.6 


7.81 


2.1 


3.3 


6.8 


4.73 


2.99 


1.45 




SW^'^ 


217 


14.9 


7.7 


28.6 


5.45 


2.3 


3.9 


6.5 


2.95 


1.78 


1.07 


MACS J0242. 5-2132 




21 


6.1 


5.5 


12.9 


1.88 


1.3 


3.0 


5.2 


1.85 


0.81 


0.46 




SW^'i 


215 


10.9 


8.0 


19.2 


7.08 


1.9 


3.3 


6.6 


4.66 


2.72 


1.37 


MACS J0429.6-0253 




40 


5.9 


5.8 


15.7 


1.19 


1.5 


3.3 


4.7 


1.02 


0.45 


0.32 






255 


6.0 


5.2 


15.5 


0.98 


1.5 


3.3 


4.6 


0.85 


0.38 


0.27 


MACS J0547.0-3904 




65 


8.5 


11.5 


15.0 


3.22 


1.7 


2.5 


5.8 


2.35 


1.63 


0.71 


MACS J0913. 7+4056 




335 


25.2 


25.6 


30.1 


120.0 


3.0 


2.2 


7.0 


50.1 


68.5 


21.6 




SE'^'^ 


155 


23.7 


17.2 


46.1 


30.0 


4.2 


5.4 


9.6 


9.01 


7.06 


3.94 


MACS J0947.2+7623 


NE^'i 


60 


20.9 


16.3 


24.0 


89.1 


2.1 


1.7 


5.1 


53.0 


68.0 


22.0 




SW^'i 


250 


18.8 


13.7 


24.6 


56.2 


2.2 


1.8 


4.8 


32.6 


39.6 


14.9 


MACS J1411.3+5212 


SE'^'^ 


135 


28.9 


22.7 


36.6 


49.1 


3.3 


1.8 


4.9 


19.1 


35.6 


12.7 


MACS J1423.8+2404 


E-'2 


90 


12.6 


16.5 


23.7 


38.7 


2.3 


4.2 


9.0 


21.3 


11.7 


5.49 






280 


14.0 


11.9 


19.3 


22.5 


1.9 


3.2 


8.5 


15.2 


8.83 


3.36 


MACS J1532.8+3021 




310 


17.4 


15.4 


37.6 


31.9 


3.5 


4.7 


8.0 


11.5 


8.59 


5.07 


MACS J1720.2+3536 


SE'^'i 


150 


4.8 


5.3 


6.7 


0.84 


0.58 


3.4 


9.4 


1.82 


0.31 


0.11 




N2 


355 


12.1 


15.0 


19.0 


16.8 


1.6 


6.1 


15.4 


13.0 


3.50 


1.39 


MACS J1931.8-2634 




90 


15.6 


12.8 


29.5 


41.7 


2.6 


1.5 


2.8 


20.7 


36.0 


19.0 






270 


15.6 


12.8 


20.6 


41.7 


1.8 


1.0 


2.8 


29.6 


51.6 


19.0 


MACS J2046.0-3430 


Sc,l 


178 


6.8 


13.7 


11.6 


7.46 


1.2 


4.1 


10.9 


7.73 


2.30 


0.87 




N2 


358 


17.4 


12.3 


23.6 


15.4 


2.4 


6.3 


16.2 


7.88 


3.75 


1.46 


MACS J2140.2-2339 


Sc,2 


180 


7.1 


6.8 


14.4 


3.43 


1.4 


2.4 


4.5 


3.06 


1.79 


0.98 


MACS JOlll. 5+0855 


N2 


350 


32.8 


33.2 


77.5 


56.5 


6.5 


7.2 


11.5 


11.0 


9.97 


6.22 




SE2 


90 


28.0 


23.7 


60.7 


24.5 


5.1 


6.1 


10.2 


6.12 


5.10 


3.05 


MACS J0257.1-2325 


NE2 


58 


17.5 


15.3 


39.5 


16.2 


2.5 


3.6 


5.9 


8.08 


5.68 


3.51 


MACS J1359. 1-1929 


N2 


345 


23.9 


22.4 


53.9 


35.3 


5.2 


16.2 


26.6 


8.56 


2.77 


1.68 




S2 


175 


11.5 


10.6 


44.1 


3.75 


4.3 


19.1 


18.3 


1.11 


0.25 


0.26 


MACS J1359.8+6231 


NE2 


70 


4.2 


3.2 


9.3 


0.14 


0.85 


3.4 


5.4 


0.20 


0.052 


0.032 


MACS J1447.4+0827 


N2 


5 


21.0 


20.4 


26.9 


123.4 


2.4 


3.5 


10.3 


64.0 


44.2 


15.1 


MACS J2135. 2-0102 


NW2 


335 


27.2 


28.2 


75.3 


42.5 


5.1 


12.9 


17.7 


10.6 


4.18 


3.04 


MACS J2245.0+2637 


NW2 


325 


8.5 


7.8 


15.5 


1.8 


1.3 


2.9 


5.9 


1.76 


0.79 


0.39 



© 2011 RAS, MNRAS 000, ??-[2T] 



6 J. Hlavacek-Larrondo, et al. 



4 CAVITY ENERGETICS 

We estimate the ene rgy stored within each of the cavi- 
ties using Eg [ 2] (e.g . iBirzan et al.ll2004l: IPunn et al.1 l2005l : 
iRaffertv" et al.ll2006l : IPunn FabianI l2006l . 120081 ). Here, p 
is the thermal pressure of the ICM at the radius of the 
bubble and estimated from X-ray data, V is the volume 
of the cavity and for a relativist ic fluid 71 = 4/3, therefore 
-£^bubbie — ^yV (thi s is also supported observationally, see 
iGraham et alll2008l l. 



^bub 



71 



71 - 1 



pV 



(2) 



We assume that the cavities have a prolate shape. The 
volume is then given hy V — 4:7rR^Ri/3, where Ri is the 
semi-major axis along the direction of the jet, and is the 
semi-major axis perpendicular to the direction of the jet. 
We define "jet" as the line that connects the nuclear X-ray 
source to the centre of the cavity. The dimensions of the 
cavities we derive are given in Table [S] 

We determine the local thermal pressure by selecting 
a set of annuli containing roughly the same signal-to-noise 
and centred on the X-ray peak. For some clusters, the data 
quality was sufficient to allow ^ 4000 counts per annuli, but 
for others we could only consider annuli with ~ 900 counts 
each. Once the regions were selected, we chose a background 
region for each cluster within the same chip but far from 
any cluster emission. Since all of our clusters are above a 
redshift of 2; > 0.3, they only occupy a small portion of 
the chip, allowing ample area for selecting a background 
within the same chip. We then proceeded in deprojecting 
the data using the Direc t Spectral Deprojection method of 
ISanders FabianI (|2007l ). For each spectrum, we fitted an 
absorbed MEKAL model using statistics and left the tem- 
perature, abundance and normalization parameter (related 
to the electron density) free to vary. Here, the absorption ac- 
counts for Galactic a bsorption and we keep it frozen at the 
iKalberla et al.l (|2005l l value. In some annuli, it was difficult 
to constrain the abundance. In this case, we kept the abun- 
dance value frozen at the cluster average which we deter- 
mined by selecting a region within the entire cluster (r < 200 
kpc). 

For MACS JOlll.5+0855 and MACS J2135.2-0102, the 
count rate was so low that we could only extract regions 
containing ^ 400 counts each, and the data quality was too 
poor to correct for deprojection effects. Here, we proceeded 
in the same way, but did not deproject the spectra and used 
C-statistics instead. For MACS 1411.3+5212 (3C 295), the 
central jet axis (r < 25 kpc) and nucleus are dominated 
by non-thermal emission. During the spectral analysis, we 
therefore remove the regions along the jet axis, so that we 
are only left with the thermal emitting gas. In all systems, we 
start the inner-most annulus at r = l'' (which corresponds 
to the Chandra Point Spread Function, PSF), so that the 
central point source is not included. 

Following this technique, we were able to derive the 
temperature (kT) and electron density (rie) as a function 
of radius for each cluster. We then calculate the deprojected 
electron pressure (pe = UekT) and entropy {S — kTrie'^^^) 
profiles, but only projected quantities for the 4 clusters men- 
tioned in the previous paragraph. Our results are shown in 
Fig. lBll of Appendix B. When deriving the cavity energetics. 



we use the local thermal pressure of the gas as an estimate 
of p in Eq.n(p = n/cT - 1.92 x UekT). 

Note that deprojection relies on the cluster having a 
"relaxed" morphology, i.e. that it can be approximated 
as being "spherical". Ten of the systems in which we find 
surface-brightness depressions are part of the X-ray bright- 
est MACS clusters and have been regarded as having a "re- 
laxed" morphology since they either have a pronounced cool 
core and very good alignment with a single central dominant 
cD galaxy or a good optical/X-ray alignment and concen- 
tric contours (lEbeling et al.1 l20ld MACS J0159.8-0849, 
MACS J0242.5-2132, MACS J0429.6-0253, MACS J0547.0- 
3904, MACS J0947.2+7623, MACS J1532.8+3021, 
MACS J1720. 2+3536, MACS J1931.8-2634, MACS J2140.2- 
2339 and MACS J2245.0+2637). Two others are part of 
the 12 most distant MACS clusters (|Ebeling et al.1 l2007l . 
MACS J0257. 1-2325 and MACS J1423.8+2404), and have 
also been considered as having a "relaxed" morphology. We 
can therefore approximate these 12 systems as being spher- 
ical and deproject the spectra. The re r aainin g 8 systems 
are not part of the lEbeling et al.l (|2007l . I2OIOI ) studies and 
have not yet been classified. MACS J2046. 0-3430 shows a 
slight offset between the central X-ray peak and galaxy, but 
the cluster still has a pronounced cool core and concentric 
X-ray contours. We therefore choose to pursue with spectral 
deprojection for this cluster. All of the remaining clusters 
have concentric X-ray contours and good optical/X-ray 
alignment, apart for MACS J2135. 2-0102. However, as 
mentioned in the previous paragraphs, the count rate was 
so low in this cluster that we could not deproject the 
spectra and used the projected quantities to calculate the 
cavity energetics. 



5 CAVITY AGES 

The power capabilities of the cavities are determined by 
dividing the energy in the bubble with its age. The 
latter is given by the buoyancy rise time, the refill 
time or the sound cr ossing time. See respectively Eq. [3] 
(IChurazov et al.ll200lh. ^a. H (McNamara et al. 2000) and 
Eq. [5] (|McNamara NulsenI [20071). Here, R is the dis- 
tance from the central BCG to the middle of the cav- 
ity (projected), S is the cross-sectional are a of the bubble 
(S = Tri^w), Cd=0.75 is the drag coefficient (jChurazov et al.l 
l200lh . g is the local gravitational acceleration {g = GM{< 
R)/R'^), r is the bubble radius (r = (i?ii?w)^^^ for an el- 
lipsoidal bubble), and Cs is the sound crossing time (cs = 
y^72/cT/ ( /im j^), where kT is the plasma temperature at 
the radius of the bubble, 72 = 5/3 and /x = 0.62). 



^bu 



R 



2gV 



trefill 



t 

L-cs — 

Cs 



(3) 



(4) 



(5) 



The buoyancy rise time is the time it takes a bubble to 
reach its terminal buoyancy velocity and depends on drag 
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Figure 2. See caption on next page. 
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MACS J2140.2 



MACS J0111.5 



Cavity 



MACS J0257.1 



I 

Potential Cavity 



MACS J1 359.8 



Potential Cavity 




10 arcsec = 45.9 kpc Potential Cavities arcsec = 60.1 kpc 




MACS J1 359.1 



Potential Cavities 



10 arcsec = 67.4 kpc 
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Figure 2. Images of the MACS clusters with cavities. There are two clusters in each horizontal line, and for each cluster we sh ow the 
0.5 — 7 keV X-ray image along with a HST image, and if not available, the R passband images from lEbeling et aTl (|2007l , 120101) taken 
with the UH 2.2-m telescope (MACS J2046. 0-3430 and MACS JOlll. 5+0855) or the Second Palomar Observatory Sky Survey image, 
POSSII, for MACS J1447. 4+0827. The 10^' scale for each object is shown with the vector in the optical images, and the black cross 
symbol shows the location of the BCG. The first 13 clusters have "clear" cavities, while the remaining 7 only have "potential" cavities 
(see "potential cavities" indicated in the figure). We emphasize that deeper observations are needed to confirm if the "potential" cavities 
are real. For the clusters with published radio maps, we plot the radio axis with a green arrow. 



forces. This is a good estimate of a bubble's age if it has 
clearly detached from its AGN and has risen, which is the 
case for most of our cavities. The refill time is the time it 
takes a bubble to rise buoyantly through its own diameter 
starting from rest. The sound crossing time is the time it 
takes a bubble to travel the projected distance from the 
central AGN to its current location at the sound speed. The 
latter is used under the assumption that the bubbles travel 
at subsonic speeds. Between the three time estimates, it is 
still not clear which is the best to use, although they should 
not vary significantly. 

The values we derive in the spectral deprojecting tech- 
nique are used to estimate the thermal properties of the 



gas (see Fig. IBll in Appendix B). To determine M(< R), 
we assume hydrostatic equilibrium (Eq. [6|, where G is the 
gravitational constant, pgas ^ 1-92 /j^m^Ue and dp/dr is the 
pressure gradient. 



M(< r) 



-1- 



G pgas dr 



(6) 



Table [3] shows the values we obtain for the power stored in 
each cavity using tcs, tbuoy and trefiii- 
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Table 4. Cooling flow properties of the MACS clusters with cav- 
ities - (1) Name; (2) Bolometric X-ray luminosity (0.01 — 50 keV) 
estimated from the Chandra data; (3) Cooling time at r = 50 kpc; 
(4) Cooling radius estimated as the radius at which tcool = 7.7 
Gyrs; (5) Bolometric X-ray luminosity (0.01 — 100 keV) of the 
clusters within the cooling radius. "^The Chandra data of these 
clusters were too poor to deproject the spectra. In this case, we 
deprojected the surface brightness proflles to obtain an estimate 
of the deprojected rie and Lx- Using these values, combined with 
the average cluster temperature within r < 200, we then deter- 
mined fcool- 



(1) 


(2) 


(3) 


(4) 


(5) 


Cluster 


Lx 


icool (50kpc) 




Lxi< J^cool) 


name 


(0.01 - 100 keV) 
1044 ergs-i 


Gyr 


kpc 


(0.01 - 100 keV) 
1044 ergs-i 


MACS J0159.8-0849 




3.0 


86 


11.8i«-j 


MACS J0242.5-2132 


21.li«-J 


2.3 


110 


16.5i«j 


MACS J0429.6-0253 




3.1 


105 




MACS J0547.0-3904 


a q+0.3 


2.7 


100 


4.8i«;^ 


MACS J0913.7+4056 


19.3+^-t 


2.1 


107 


15.2i«-J 


MACS J0947.2+7623 


37.6i«-> 


1.6 


120 


29 Q+O-6 
^^•^-0.6 


MACS J1411. 3+5212 




4.8 


60 


4.oi«:^ 


MACS J1423.8+2404 


23.3+«-j 


2.1 


95 


16.3i«-J 


MACS J1532.8+3021 


33.2+-^-^ 


1.5 


115 


24.9t«-J 


MACS J1720.2+3536 


15.li«;l 


3.8 


100 




MACS J1931. 8-2634 


9.0+0-} 


1.8 


112 




MACS J2046.0-3430 


13.4+«-4 


3.1 


81 


o O+0.3 
°-"^-0.3 


MACS J2140.2-2339 




1.9 


107 


14 2+^-^ 


*MACS JOlll. 5+0855 


^ O + 0.6 


5.7 


78 


9 O+0.5 
^••^-0.3 


MACS J0257. 1-2325 


17A+1-1 


8.0 


54 


o 4+0.6 


MACS J1359.1-1929 


' -^-0.3 


3.9 


76 


4.2l«;2, 


MACS J1359.8+6231 


a O + 0.4 

°-°-0.3 


5.9 


67 


2.9l«;} 


MACS J1447.4+0827 


44.0 tli 


1.3 


125 




*MACS J2135.2-0102 




11.9 






MACS J2245.0+2637 


11 g+0.6 


4.0 


89 


6 3+^-4 



6 COOL CORES 

6.1 Clusters in the MACS sample with cavities 

Statistical studies of cavities in clusters of galaxies have 
shown that cavities sit predominately in cool core clusters 
(jPunn et al.ll2005l : IPunn Fabianll2006h . In order to com- 
pare our results with theirs, we have computed detailed cool- 
ing time profiles for all of our clusters with cavities. More 
precisely, we compute the cooling time (tcooi) by using Eq.[7l 
and the thermal gas properties derived in Section 4 (shown 
in Fig.lBB. 

5 1.92 ue kT V 

tcool = (7j 

Here Ue is the electron density, kT is the gas temperature, 
Lx is the gas X-ray luminosity and V is the gas volume con- 
tained within each annulus. We use the deprojected spectral 
quantities as an estimate of the plasma parameters and show 
our results in Fig. [S] However, we show in black the two 
clusters for which we were not able to deproject the spec- 
tra (MACS JOlll.5+0855 and MACS J2135.2-0102). In this 
case, we deprojected the surface brightness distribution, and 
then converted the counts per second of each annulus into 
the predicted electron density and flux of a mekal model 
for which the temperature and abundance were frozen at 
the average value within r < 200 kpc. We then used these 
values of Ue and Lx combined with the average cluster tem- 
perature within r < 200 kpc to determine tcooi using Eq. 
[71 This is possible since X-ray observations in clusters give 
very precise estimates of the surface brightness profiles. It 
is therefore much easier to deproject the surface brightness 



profiles than the spectra, allowing us to obtain an estimate 
of the "deprojected" cooling time profiles for these two ob- 
jects. 

Fig. [3] shows that many of the clusters in which we 
detected cavities are strong cool cores, with cooling times 
reaching down to ^ 0.5 Gyr in the central regions. The esti- 
mates we obtain for the central cooling times are shown in 
Table [H We define the central cooling time as the cooling 
time we derive at a radius of 50 kpc and not the cooling 
time of the inner most bin. This is because the radius of the 
inner most bin depends not only on the cluster distance, but 
more importantly on the data quality, and is therefore not 
neces sarily the b est estimate of the central cooling time (see 
also Bauer et al. 1[2005). For almost all of the clusters in our 
sample we have cooling time estimates below 50 kpc, allow- 
ing us to obtain an accurate central cooling time at 50 kpc. 
For the remaining clusters, we determine a rough estimate 
by extrapolating the values down to 50 kpc. 

Adopting this definition of the central cooling time, we 
find that 19 of the 20 clusters in which we identified cavi- 
ties are mild cool cores with tcooi < 8 Gyrs and that all of 
the systems with "clear" cavities are strong cool cores with 
tcooi < 3 — 5 Gyrs. We also illustrate this in the top panel 
of Fig. m where we plot the distribution of cooling times in 
our sample of clusters with cavities. 

Since many clusters dX z — 1 appear to have similar 
properties to present-day ones, they can be regarded as re- 
laxed systems and more importantly, we can assume that 
a cooling flow should have had the time to establish itself 
since then. We therefore choose to define the cooling radius 
of each cluster as the radius at which the cooling time is 
equal to the z = 1 look-back time. For our cosmology, this 
corresponds to 7.7 Gyrs. This definition of th e cooling ra- 
dius was also adopted bv lRaffertv et al. liooi), 

allowine: us 

to directly compare our results with theirs. Although there 
are oth er statistical studies of cavitie s in the literature (e.^ . 
Birzan et al. 2004; Dunn et al. 2005; Dunn Fabianll2006l : 
Birzan et a l. 2008; Dunn Fabian 2008h, many of the ob- 
jects have already been included in iRaffertv et al.l (|2006h . 

Finally, we derive the bolometric X-ray luminosity of 
each cluster (Lx (0.01-100 keV)), and the cooling luminos- 
ity which we define as the bolometric X-ray luminosity 
(0.01 — 100 keV) within the cooling radius, Lx(< Tcooi). 
When deriving these quantities, we exclude the central point 
source (r < 1^^) for clusters with point-like nuclei at their 
centres, since some have important non- thermal emission 
associated with them. These include MACS J0547.0- 
3904, MACS J0913.7+4056, MACS J0947.2+7623, 
MACS J1411.3+5212, MACS J1423.8+2404, 
MACS J1931.8-2634 and MACS J2046.0-3430. 



6.2 MACS cluster sample 

In order to derive the fraction of cool core clusters with cav- 
ities in the MACS sample, we computed the central cooling 
time within 50 kpc for all 76 MACS clusters with Chandra 
archival data. Here, we extracted a spectrum for each cluster 
within r < 50 kpc and fitted a mekal model to the spectrum 
while letting the temperature, abundance and normalisation 
parameter free to vary. For some clusters, the counts were 
not sufficient to constrain the fit. In this case, we kept the 



© 2011 RAS, MNRAS 000, ??-[2T] 



10 J. Hlavacek-Larrondo, et al. 



abundance and temperature frozen at the value derived for 
the entire cluster (r < 200 kpc). 

We could not compute detailed cooling time profiles for 
all the MACS clusters since the data quality in many is not 
sufficient to derive deprojected or even projected profiles. 
We therefore choose to use the cooling time within the cen- 
tral 50 kpc (tcooi(^ < 50kpc)) as an estimate of the central 
cooling time. This definition of the central cooling time is 
not the same as the one used previously and defined as the 
cooling time derived at a radius of 50 kpc (tcooi(r = 50kpc)). 
To determine the fraction of cool core clusters with cavities 
in the MACS sample we therefore recalculated the central 
cooling time of clusters with cavities adopting the new defi- 
nition (tcooi(^ < 50kpc)). Our results are shown in the bot- 
tom panel of Fig. [H 

We find that 37 of the 76 MACS clusters have a strong 
central cool core (tcooi(^ < 50kpc) < 3 Gyrs ), corresponding 
to a cool core fraction of 50 per cent. Bau er et al.l (|2005l ) 
derive detailed cooling time profiles for 38 X-ray luminous 
clusters from the ROSAT BCS survey (0.15 <z< 0.4) and 
find that 11 have a central cooling time less than 3 Gyrs 
(and 55 per cent less than 10 Gyrs). They use the cooling 
time derived at a radius of 50 kpc as their definition of the 
central cooling time, whereas we use the cooling time de- 
rived within 50 kpc for the central cooling time which gives 
smaller cooling times. Nonetheless, we find that a significant 
fraction in our sample are cool core clusters and out of these 
37 clusters, we find that 19 have cavities, corresponding to a 
cavity detection rate in cooling systems of 50 per cent. Note 
however, that more than half of the remaining cool cores 
(^cooi(r<50kpc) < 3 Gyrs) with no cavities have an exposure 
time below 30 ks (see top panel of Fig. [5]), which could ex- 
plain why we do not detect cavities in these systems. We 
also plot the central cooling time as a function of Chandra 
exposure time for all MACS clusters in the bottom panel of 
Fig. O and discuss these results in Section 9.1. 



7 CAVITY SCALING RELATIONS 

We now address the various scaling relations known for the 
cavity population in groups and clusters. It has been known 
for quite some time that AGN feedback is finely tuned to 
the energy needed to offset cooling of the ICM. If the input 
of energy were too low, then the core would cool rapidly and 
cause massive quantities of gas to condense and form stars 
at rates much higher than what we observe. However, if the 
input of energy were too high, then this would cause the 
core to over-heat and we would not observe the low central 
cooling times in the cores of some clusters. 

To illustrate the fine tuning of AGN feedback, we plot 
the mechanical power arising from the central AGN as a 
function of the bolometric X-ray luminosity within the cool- 
ing radius (Lcooi(0.01 - lOOkeV)) in Fig. [1 Within this ra- 
dius, the gas should have had the time to cool if no addi- 
tional heating sources were present. Although there is some 
cooling happening in the cores of clust ers, studies suggest 
it is only on the or der of 10 per cent (|Birzan et al ] l2004l : 
iRaffertv et al]l2006l ) or even less than 10 per cent accord- 
ing to th e more recent stud i es wi th high-resolution spec- 
troscopy ([Sanders et al ] |2008l . l2010h . and therefore does not 
affect significantly the relation plotted in Fig. [G] The me- 




10 100 
Radius (kpc) 

Figure 3. Cooling time profiles for all clusters in our sample for 
which we were able to identify cavities. Each colour illustrates 
a different cluster. The clusters shown with the black lines are 
those in which the Chandra data were too poor to deproject the 
spectra. In this case, we deprojected the surface brightness profiles 
to obtain an estimate of the deprojected rie and Lx- Using these 
values, combined with the average cluster temperature within r < 
200, we then determined tcool 

. MACS J2135. 2-0102 is the only object in our sample with 
cavities that has a central cooling time larger than 10 Gyrs 
(shown with the black dotted line). 

chanical power is estimated using the energy stored within 
the cavities of each system {4pV) and the buoyancy rise 
time (tbuoy)- We also include the cavities of IRaffertv et al.] 
(|2006l ) and illustrate in red the two systems common to both 
samples (MACS J1423.8+2404 and MACS J0947.2+7623). 
For MAC S J1423 .8+2404, our results agree weh with those 
of Raffertv et al ] ([2OO6I but there is a large discrepancy 
be tween our estimate o f the mechanical power and that 
of IRaffertv et al.1 (|2006l ) in MACS J0947.2+7623. How- 
ever, the more r ecent study on MACS J0947.2+7623 by 
Cava gnolo et al.l (| 20 111 ) finds an outburst power consis- 
tent with our estimate of Lmech (^ 6 x 10^^ erg s~"^ ). We 
do not include the data po ints of iDunn et al l (|2005h or 
iDunn FabianI (|2006l . [2008') as they define the coohng ra- 
dius as t he radius at wh ich tcooi = 3 Gyrs, nor the data 
points of Cava gnolo et al.l (|2010( 1 as they only calculate the 
cooling luminosity within 0.3 — 2.0 keV. Fig. \6\ shows that 
the MACS sample contains some of the most powerful out- 
bursts known, and that even at these extreme energies, the 
central AGN is able to offset cooling of the ICM at a fine 
level. 

We also look at the properties of the mechanical power 
arising from the central AGN as a function of radio luminos- 
ity. As these cavities are being created by the central AGN 
and being filled with relativist ic electrons, one would expect 
a correlation between the mechanical power and radio lumi- 
nosity, which is a measure of the radio power arising from 
the AGN. This has been looked at in detail for clusters, el- 
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Figure 4. Top panel: Distribution of ^cool cit 50 kpc for the MACS 
clusters in which we identified cavities. Here, we derived the cen- 
tral cooling time by extracting a series of annuli and determining 
the cooling time at r = 50 kpc. Bottom panel: Distribution of 
tcool within 50 kpc (r < 50 kpc) for all 76 MACS clusters with 
Chandra observations. Among these, we outline those with clear 
bubbles in black (all of these have a radio source associated with 
the central galaxy) and those with potential cavities and a radio 
source in grey, as well as those with potential cavities but no radio 
source in diagonal lines. 



lipticals an d more recently groups of gjalaxies (iBirzan et alJ 
l2004l . l2008l : ICavagnolo etaLl I2OIOI : lO ' Sullivan et al.ll201lh . 
although the scatter remains large and extends over several 
orders of magnitude (> 3). We have computed the radio lu- 
minosities for our sample (see Table[2]) and show in Fig. [Tithe 
mechanical power arising from the central AGN {4pV/thuoy) 
as a function of radio luminos i ty (L radio)- We include the 
data points of 'O' Sullivan et al.' ("20111) for galaxy groups, as 
well as those in Rafferty et al. (2006[j_, but use the estimate 
of Lradio given in Birzai^t al.l (I20p8 | l. W e do not include 
the data points of ICavagnolo et al.l |2010l l as they do not 
include an estimate of Lradio in their paper. We discuss this 
plot in Section 9.2. 



8 RESULTS 

In this section we review the results of each system. We 
begin by addressing the systems with "clear" cavities, and 
then address those with "potential" cavities. Many of the 
MACS clusters have been studied in terms of their mass 
distribution and cosmological implications (e.g. I Allen et"aI1 
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Figure 5. Chandra exposure time distribution for the 37 MACS 
clusters with tcool(r<50kpc) < 3 Gyrs (top panel) and for the 

> 3 Gyrs (middle panel). 



39 MACS clusters with t 



cool(r<50kpc) 



Among these, we outline the 13 with clear bubbles in black (all of 
these have a radio source associated with the central galaxy) and 
those with potential cavities and a radio source in grey, as well 
as those with potential cavities but no radio source in diagonal 
lines. Note that more than half of the cool cores (^cool(r<50kpc) 
< 3 Gyrs) with no cavities have an exposure time below 30 ks, 
which could explain why we do not detect cavities in these sys- 
tems. Also shown in the bottom panel is the central cooling time 
as a function of Chandra exposure time. The MACS clusters with 
clear bubbles are shown in dark blue, and those with potential 
bubbles are shown in light blue. The straight line shows the cut- 
off between those with tcool(r<50kpc) < 3 Gyrs (light grey circles) 
and those with ^cool(r<50kpc) > 3 Gyrs (black circles). 
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Figure 6. Mechanical energy being injected into the ICM by the 
central AGN (estimated from the energetics of the cavities) as 
a function of energy needed to offset cooling of the ICM by the 
central galaxy. Those common to both samples are shown in red 
and yellow. 
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Figure 7. Mechanical energy being injected into the ICM by the 
central AGN (estimated from the energetics of the cavities) as a 
function of total integrated radio luminosity (lOMHz to lOGHz) 
of the central source. Those common to both samples are shown 
in red and yellow. 



I2OO8I : lEttori et all l2009l : IPuetzfeld l2005l ). as well as 

in terms of h ow the cluster met a llicity content evolve s 
with redshift (jBalestra et al ] L2OO7I; i Maughan etHI 120081 ). 
These include MACS J0159.8-0849, MACS J0242.5- 
2132, MACS J0429.6-0253, MACS J0947.2+7623, 
MACS J1423.8+2404, MACS J1532.8+3021, 
MACS J1720.2+3536, MACS J1931.8-2634, MACS J2046.0- 
3430, MACS J2140. 2-2339 and MACS J1359.8+6231. 

8.1 Clusters with "clear" cavities 

8.1.1 MACS J0159. 8-0849 

We find two clear X-ray cavities in this system. They are 
located symmetrically opposite with respect to the central 
AGN, and at a distance of ~ 25 kpc from the core. The 
power stored within them (Pbuoy ^ 5 x 10^^ erg s~"^ ) falls 
short by a factor of 2 from the cooling luminosity. The HST 
optical image shows interesting features, including a poten- 
tial double or triple nucleus with streams of material leading 
back to the central nucleus. 



8.1.2 MACS JO242.5-2132 

We find two cavities in this cluster, with the south-west 
one being clearer and somewhat larger than its northern 
counterpart. The optical image shows that there is a bright 
central isolated galaxy at the centre, which coincides with a 
very strong radio source of flat spectrum (~ 1 Jy bright, the 
second strongest in our sample). The power stored within 
the cavities (Pbuoy ^ 4 x 10^^ erg s~^ ) falls short by a factor 
of 4 from the cooling luminosity. This cluster is part o f the 
34 most luminous MACS clusters (|Ebeling et al.ll2010l ). 



8.1.3 MACS JO429. 6-0253 

This cluster has two clear cavities being inflated by the cen- 
tral AGN. The cavities are quite small (r ^ 5 kpc), and 
the power stored within them falls short by almost an or- 
der of magnitude when compared to the cooling luminos- 
ity (Pbuoy ^ 10^^ erg s~^ ). There is a hint of another po- 
tential cavity to the north-east, but deeper X-ray observa- 
tions are needed to confirm this. The optical image shows a 
dominant central galaxy harbouring a bright nucleus. This 
cluster is also part o f the 34 most luminous MACS clusters 
(jEbeling et al.ll2010l ). 

8.1.4 MACS JO54 7. 0-3904 

We find one cavity in this cluster of radial dimensions 
8.5 X 11.5 kpc^, surrounded by bright rims. The power stored 
within it alone (Pbuoy 2 x 10^^ erg s~^ ) can prevent the 
gas from cooling to almost half of the cooling luminosity. The 
dominant central galaxy harbours a bright point-like nucleus 
which shows hints of a jet in the optical, aligned with the jet 
axis along the X-ray cavity. The X-ray image also shows a 
bright nucleus dominating in emission. This cluster is part o f 
the 34 most luminous MACS clusters (jEbeling et ahlboiOl ). 

8.1.5 MACS J0913. 7+4056 

The central galaxy in this cluster, more commonly known 
in the literature as IRAS 09104+4109, is a hyper-luminous 
infrared cD galaxy (Lir ^ 12 x 10^^ L©; iKleinmann et all 
1988) hosting a heavily ob scured Type 2 AGN (- l-5x 10^^ 
cm "^: IVignali et al.ll201lh . The Chandra X-ray spectrum of 
the nucleus shows that there is a prominent Fe Ka lin e at 
the redshift of the source (see also llwasawa et al]|200lh . 
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We find two cavities in this system. The cavity t o the 
north-west was first noticed by llwasawa et alj (|200lh , but 
the 9.1 ks Chandra observations were too faint to confirm 
whether this cavity was real. By using a much deeper ob- 
servation of 70.4 ks, we confirm the existence of the first 
cavity and find evidence of a second located at a simi- 
lar radius to the south-east. Although larger, the cavities 
carved out in the X-ray image coinci de with the location 
of the radio hot spots seen at 20 cm (IHines Willslll99l 
also see Fig. [2] where we show the radio jet axis), and the 
optical image shows a spectacular set of filaments spread- 
ing towards the north-east and south. These filaments run 
along both cavities that we have identified. However, one of 
these filaments extends further out (rmax(fiiament) = 51 kpc 
and rmax(cavity) — 55 kpc) and is directed more eastward 
than the northern cavity. As seen in many nearby clusters , 
including M87 in the Virgo cluster teohringer et al.l ll995l : 
IChurazov et~alll200lh and NGC 1275 in the Perseus Cluster 
(|Fabian et al.ll2003l \ filaments usually trail behind cavities, 
and could indicate that there is another, older cavity to the 
north-east, although it is difficult to tell with the current 
X-ray observations. The power stored within the cavities is 
substantial (Pbuoy 8x 10^^ erg s~^ ) and is more than suffi- 
cient to prevent the ICM from cooling. This cluster harbours 
the third most powerful outburst in our sample (in terms of 
total power), and has the cavity with the largest enthalpy 
in those classified as "clear" cavities (pV ^ 10^*^ ergs~^). 
MACS J0913.7+4056 has one of the most powerful AGN 
outbursts known in clusters. 



8.1.6 MACS J 094 2+ 7623 

MACS J0947.2+7623 (RBS 0797) is another cluster known 
to have one of the most powerful AGN outbursts (see 
ICavagnolo et al.l 2011, for a detailed analysis on the AGN 
outburst energetics). It has a very strong cool core, re- 
quiring extreme mechanical feedback from its central AGN 
(Lcooi 3 X 10^^ erg s~"^ ), and the X-ray image shows 
two large X-ray cavities surrounded by bright rims, as well 
as a point-like nucleus. This cl uster is also part o f the 34 
most luminous MACS clusters (jEbeling et ahlboiol v In the 
optical, we see a large central galaxy hosting a complex 
and disturbed core. The power stored within each cavity 
(i^buoy 4 — 7 X 10^^ erg s~"^ ) is more than sufficient to 
prevent the surrounding gas from cooling. This cluster has 
the most powerful outburst in our sample and the e nerget - 
ics we derive are consistent with Cavagnolo et al.l (|201lh . 
except that we estimate the size of the north-east cavity to 

be slightly larger.^ 

Interestinglv. lGitti et aP (|2006l ) found that the 4.8 GHz 
radio jet axis is perpendicular to the jet axis running through 
the cavities (see Fig. [2] for the direction of t he 4.8 GHz radi o 
jet a xis). This was later confirmed bv Birzan et al.l (|2008l V 
I Cavagnolo et al.l (|201lh also suggested that there is a hint 
of a small, newly formed cavity at the tip of the southern 
4.8 GHz jet. These authors furthermore find that the bright 
rims surrounding the cavities could be driving weak shocks 
into the ICM and suggest that the unusually deep surface 
brightness depressions of the cavities, combined with the 
ambiguous correlation between the X-ray and radio mor- 
phologies, is more consistent with these cavities being elon- 



gated along the line of sight. These cavities could therefore 
be larger than measured in the plane of the sky. 

8.1.7 MACS J1411-3+5212 

MACS J1411.3+5212 (3C295) hosts a very powerful cen- 
tral AGN embedded in a giant cD galaxy. This galaxy 
has a bright radio core and is one of the most power- 
ful FRII sources known. The cluster shows evidence of 
non-thermal emission filling the in ner radio lobes , con- 
sistent with ICCMB photons (e.g . iBrunetti" et al.1 I2OOII : 
lAllen et al.ll200ll : iHarris et al.ll200Ql i. The radio maps show 
GHz radio emission filling the inner ~ 25 kpc of the 
core, with a j et axis aligned in the north-west to south- 
east direction (|Perlev Tavloilll9"9ll:lTaylor Perley|[l99i ; 

iHarris et alTl200C 



IBrunetti etHI I2OOII : fAllen et al.1 12OOII : IHarris et al.1 lioool . 
see also Fig. [2 | for th e direction of the GHz radio jet axis). 
IBrunetti et al.l (j200l|) also find a prominent Fe K line at 
the redshift of the source associated with the nucleus of the 
central galaxy, consistent with this object being a powerful 
hidden and absorbed AGN. 

During the spectral analysis of this source (Sections 4, 
5 and 6), we omitted the central regions containing the IC- 
CMB X-ray emission since our aim was to extract the prop- 
erties of the thermal gas. The cavity that we find is located 
at the edge of south-east ICCMB X-ray lobe, and could 
be an older cavity carved in the past by the radio source 
but which has now risen buoyantly. The radio maps at GHz 
frequencies s how that the emission is limited to the inner 



r < 25 kpc (IPerlev k Taylor "l99l"; "Taylor k Perlev 'l992l: 
IBrunetti et"al] I2OOII : iXllen et al. 2001; Harris et al. 20q3). 



The cavity we find is beyond this radius (rcavity 35 kpc), 
and is consistent with being older since the GHz radio obser- 
vations do not seem to show emission extending out to this 
radius. This cavity is large and the power stored within it 
(Pbuoy ^ 4 X 10^^ erg s~^ ) is more than sufficient to prevent 
the gas from cooling. 

8.1.8 MACS JI423. 8+2404 

This cluster hosts a very large central galaxy harbouring a 
bright nucleus seen both at optical and X-ray wavelengths. 
In the optical, the HST image shows that the central galaxy 
is heavily disturbed: the bright point-like nucleus is sur- 
rounded by a thick envelope in the inner regions (r < 2.6 
kpc) and a more diffuse and disrupted envelope out to r 20 
kpc. The optical image also shows two depressions, most 
likely due to dust: one to the west which coincides with the 
location of the cavity seen at X-ray wavelengths and a hint 
of a second to the north-east of the nucleus which coincides 
with a very small depression seen in the X-rays (although we 
have not included this inner smaller cavity in our sample). 
The power stored in the cavities (Pbuoy 2 x lO'^^ erg s~"^ ) 
is sufficient to prevent the gas from cooling within the 
cooling radius. This cluster is also the most distant object 
in the MACS sample in which we have identified cavities 
(z = 0.5449), and is part of the 12 most d istant clusters 
detected by the MACS survey at z > 0.5 (Ebeling e t al.l 
2007). The bright nucleus does not show any evidence of a 
hidden quasar in the form of Fe K lines at the redshift of the 
source, as opposed to other MACS clusters mentioned previ- 
ously with a bright X-ray nucleus. iKartaltepe et al.l (|2008l l 



© 2011 RAS, MNRAS 000, ??-[2T] 



14 J. Hlavacek-Larrondo, et al. 



looked at the large scale structure of the 12 most distant 
clusters detected by the MACS survey at z > 0.5 by tracing 
the surface brightness density of galaxies near the red se- 
quence, and singled out that MACS J1423. 8+24Q4 was the 
most relaxed cluster in the sample. Limo usin et al 
also confirms that the cluster seems to be nearly fully viri- 
alized as suggested by gravitational lensing analysis. Having 
such a virialized cluster only ^ 7 Gyrs after the Big Bang 
puts constraints on structure formation and evolution in a 
cosmological context. 



8. 1.9 MACS J 1532.8+3021 

MACS J1532.8+3021 (RX J1532.9+ 3Q21) has the fourth 
most infrared luminous BCG in the iQuillen et al. 1 (|20Q8 ') 
sample of 62 X-ray luminous clusters, enough that it can 
be classified as a luminous infrared galaxy (LIRG) with 
LiK = 22.6 X 10^^ erg s~^ . The central galaxy also has 
one of the most massive molecular gas detections known 
(2.5 X 10^^ M©, the most massive detection in Edge i2001|) 
and is the second most opti cally line-lumino us central clus- 
ter galaxy in the sample of I Crawford et al with an 
optical line emission of 4.2 x 10^^ erg s~^ . The galaxy is 
also particularly blue, suggesting star formation, and has 
the largest star formation region in the lO'Dea et al.l (|2008l l 
cluster sample, with a scale length of ~ 70 kpc. The cluster 
is also particularly massive M> lO^^Mo, as found by the 
weak lensing study of IPahle et al.l (|2002 ). 

We find one clear cavity associated with the central 
galaxy, located to the west. It is quite large (r ^ 20 kpc) and 
situated almost 40 kpc from the nucleus. This cluster has 
the third strongest cooling luminosity in our MACS sample 
with cavities, requiring around 3 x 10^^ erg s~^ to prevent 
the gas from cooling. We find that the power stored within 
the cavity falls short by a factor of 2 — 3 of the cooling 
luminosity. The HST image suggests that there is some fila- 
mentary structure surrounding the central galaxy (although 
it is not as clear as in MACS0913. 7+4056), and one of the 
northern filaments seems to sweep around the bottom por- 
tion of the cavity. 



8.1.10 MACS J 1720.2+3536 

We find two cavities associated with the central regions of 
MACS J1720. 2+3536 (Z8201), one to the north and another, 
though very small and not as clear, to the south-east. If 
we consider the power stored within the cavities using the 
buoyancy time-scale, then the power falls short by a factor 
of 2 of the cooling luminosity. If we consider the sound- 
crossing time, then the power is more than sufficient to off- 
set cooling of the ICM. The HST image shows that the 
central galax y has a bright c ore su rrounded by a large dif- 
fuse envelope. iRaffertv et gH (|2008l ) also find that the central 
galaxy has blue emission extending out to 2.5^^ 14 kpc). 
This cluster is part of the 34 most luminous MACS clusters 
(jEbeling et al.ll201Ql ). 



8.1.11 MACS J 1 931 . 8-2634 

MACS J1931. 8-2634 is a strong cool core cluster harbouring 
one of the most powerful AGN outbursts known. The cen- 



tral nucleus has carved two large X-ray cavities to the wes t 
and east, greatly disrupting the core (|Ehlert et al.ll201lh . 
Diffuse 1.4 GHz radio emission fills the central regions and 
extends in the west-east direction, but no clear radio lobes 
or jets are seen. The X-ray image (Fig. [2]) shows two bright 
ridges formed to the north and south, some 25 kpc from 
the nucleus. The north ridge is also mad e of cool, metal-rich 
gas rich in Ha filaments and young stars. lEhlert et all (|201lh 
propose that these ridges were formed by the AGN outburst, 
as well as sloshing of the core al ong the north to so uth di- 
rection induced by a past merger. lEhlert et al ] (|201lh found 
that the temperature profile decreased with radius, but that 
in the inner-most regions (r < 30 kpc), the temperature 
increased once more. This could indicate shock heating. In 
Fig. EH we show the spectral profiles of MACS J1931.8 we 
obtained during our analysis. Our temperature profile also 
increases in the inner- most regions (r < 30 kpc). 

The cavity edges are not well defined, and the radio 
emission shows no clear radio lobes. Determining a pre- 
cise measurement of the cavity sizes is therefore difficult. 
Taking this into account, we estimate that the 4pV energy 
stored within the cavities is roue: hly - 3 X 10^° erg and the 
power input is Pbuoy 8 — 9 x lO'^^ erg s~^ or using the 
time-c rossing time, Pcs 4 — 5 x 10^^ erg s~^ . lEhlert et al.l 
(|201lh find a total energy of - 1 - 8 x 10^° erg {4pV) and 
Pcs ~ 4 — 14 X lO'^^ erg s~^ , both consistent with our results. 
MACS J1931.8-2634 therefore has one of the most powerful 
AGN outbursts known. This cl uster is also part o f the 34 
most luminous MACS clusters (|Ebeling et al.ll201Ql V 



8.1.12 MACS J2046. 0-3430 

We find two cavities associated with this cluster. The first 
is located to the south, and has a fiat shape with bright 
rims surrounding it. The second cavity is located to the 
north but is less well defined. We were not able to find 
any HST archival data and show in Fig. [2] the R pass- 
band UH 2.2-m image. The optical image shows that the 
central dominant galaxy is offset from the X-ray point 
source by ^ 2.7^^ (15 kpc). The power stored in the cavi- 
ties (Pbuoy ^ 6 X 10^^ erg s~^ ) is comparable to the amount 
of energy needed to prevent the gas from cooling. 



8.1.13 MACS J 21 40. 2-2339 

MACS J2140.2-2339 (MS2137.3-2353) is a luminous cluster, 
dominated by a large central galaxy. The most interesting 
feature seen in the optical image is a large radial lensing arc 
some 70 kpc from the nu cleus, which was the first of its kind 
to be discovered (jFort et al. 1992) and has been looked at 
in de t ail by many authors (e .g.lGavazzi et al. | 1 20031: ICavazzj 
20051 : IComerford et al.ll2006l : ISand et al.ll2008l : iMerten et al.1 
2009l l. T his cluster is part o f the 34 most luminous MACS 
clusters ('Ebeling e t "all 1 20 id ). We find a small cavity to the 
south of the nucleus, only capable of preventing the gas from 
cooling to 10 per cent of the cooling luminosity. Although 
we classify this cavity as being "clear" , it is very small and 
does not have any bright rims surrounding it. 
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Figure 8. From left to right and top to bottom: Distribution of the mechanical power being injected by the central AGN, the total 
energy output stored in the cavities, the ratio of the distance from the core and the average radius of the bubble, the average radius, the 
distance from the core and the buoyancy rise time. 



8.2 Clusters with "potential" cavities 

We now address the systems with "potential" cavities. Al- 
though we found evidence for cavities in these clusters, many 
of the X-ray observations have a low number of counts, and 
deeper observations are needed to confirm if these cavities 
are real. 

8.2.1 MACS JOlll. 5+0855 

We find two cavities associated with the dominant central 
galaxy in this cluster. They are both quite large (30 kpc 
in length) and located some 60-70 kpc from the nucleus. 
However, the data quality for this system is quite poor, and 
we were not able to deproject the spectra. Furthermore, we 
were only able to extract two annuli to derive the varia- 
tions of the thermal properties with radius (see Fig. IB1|) . 
Although the central cooling time we derive at r = 50 kpc 
is below 10 Gyrs, it is not clear whether this cluster has a 
cool core because the data quality does not allow us to ob- 
tain a detailed temperature profile in the inner regions. This 
cluster also has no central radio source associated with it. 
If we assume that this cluster has a small cool core, then 
the cooling luminosity within the cooling radius (rcooi = 78 
kpc) is ^ 2 X 10^^ erg s~^ . The cavity power is on the order 
of Pbuoy 15 X lO'^'^ erg s~"^ , but since these cavities are not 
very well-defined (and only classified as potential cavities), 
this value could be quite uncertain. 

8.2.2 MACS J0257. 1-2325 

This cluster is located at a redshift oi z — 0.5039 and 
is part of the 12 most distant clusters detected by the 



M ACS survey at z > 0.5 (jEbeling et al.ll2007l ). The study 
bv iKartaltepe et aP (|2008l l looked at the large scale struc- 
ture in the 12 most distant MACS clusters (including 
MACS J1423.8+2404), and found that MACS J0257.1-2325 
has the most complex system in terms of large-scale struc- 
ture (on Mpc scales), with infalling systems along cosmic 
filaments. The HST image in Fig. [2] shows that the clus- 
ter is dominated by a large central galaxy, hosting a bright 
core and diffuse envelope. We also find many lensed arcs sur- 
rounding the central regions (r < 70 kpc), including a very 
bright one to the south-west. We find no radio source asso- 
ciated with the central galaxy, and only a mild cooling flow 
with a central cooling time of 8 Gyrs. Although, we were 
able to deproject the spectra, it is not clear if the temper- 
ature profile is decreasing in the inner regions because the 
error bars are still quite large. The cavity we find is located 
to the north-east some 40 kpc from the nucleus and the 
power stored within (Pbuoy 6 x 10^^ erg s~^ ) is sufficient 
to prevent the gas from cooling. 

8.2.3 MACS J1359. 1-1929 

This cluster hosts a large central galaxy, surrounded by a 
diffuse envelope as seen in the HST image, and has a ra- 
dio source associated with the central regions. We also find 
a central cooling time of 4 Gyrs and a temperature profile 
steadily declining in the inner regions, with a cooling lu- 
minosity on the order of 4 x lO'^'^ erg s~^ . We were able to 
identify two cavities in this system, although not well de- 
fined and deeper data is needed to confirm their existence. 
If these cavities are real, then the power stored within them 
is Pbuoy ~ 3 X 10^^ erg s~^ . The X-ray image also shows 3 
small depressions (^ 6 kpc in size) immediately surrounding 
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Figure 9. Left - Energy of each cavity as a function its average radius (/^average)- The dotted line shows the slope of 3 one would 
expect between pVtot and /^average since pVtot cx: i^average- Middle - Energy of each cavity as a function of the pseudo-pressure, the latter 
determined by dividing pVtot by the average radius cubed (i^average^)- Right - Average radius cubed (i^average ^) of each cavi t y as a 
function of the pseudo-pressure. The light blue points are for the MACS clusters and the black points are for the iRaffertv etHI ((2OO6I ) 
objects. 



the central X-ray peak, which could also contribute to the 
heating. 



8.24 MACS J1359. 8+6231 

MACS J1359.8+6231 (MS 1358.4+6245 or EMSS 
1358+6245) is a rela xed cool core cluste r known for 
itslensinff proper t ies (iHoekstra et aL 

I Il998l : I Franx et al. 

I1997I : lAllenI Il998l: IZitrin et a" l201lh and detailed mass 
estimates (|Arabadus et al.|[2002 V 

We have identified one small cavity to the north-east 
located near the central galaxy. The optical image shows 
that the central regions are dominated by a large galaxy with 
a bright core and diffuse envelope. The power capabilities of 
the cavity are quite small with Pbuoy ^ 5 x 10^^ erg s~^ and 
are not capable of preventing the cool core from cooling by 
more than an order of magnitude. This cluster stands out in 
the Lmech — i^cooi plot showu lu Fig. [6]since the mechanical 
energy within the cavity fails to meet the power needed to 
offset cooling. 



8.2.5 MACS JI447.4+O827 

MACS J 1447.4+0827 (RBS 1429) has been recently- added 
to the MACS cluster sample. It is an ultra-luminous clus- 
ter (the most X-ray luminous in our sample with cavities; 
Lx ~ 5 X 10^^ erg s~^ ) and has a very strong cool core, 
therefore requiring extreme mechanical feedback from its 
central AGN (Lcooi -3-4x10^^ erg s"^ ). With just - 10 ks 
on the source, we find evidence of a large cavity, although 
not well defined. The energetics of the cavity in itself are 
substantial, and the power stored within this outburst is 
enough to prevent the gas from cooling. This cluster has the 
most energetic outburst in our sample in terms of enthalpy 
(^ 10^° erg), making it one of the most powerful outbursts 
known. 



8.2.6 MACS J2135. 2-0102 

MACS J2135.2-0102 (IRXS J213515.7-010208) is a lumi- 
nous cluster known for its lensing properties. An earlier 
study found a series of arcs in the shape of a "Cosmic 



Eye" and "Cosmi c Eyelash" associated with this cluster (se e 
ISmail et al.]l2007l : ISiana et al.ll2009l : ISwinbank et al.ll2010l ). 

However, the data quality of the Chandra X-ray image is 
very poor, and we were not able to deproject the spectra. We 
find no radio source associated with the central regions and 
that the central cooling time is more than 10 Gyrs. There is 
however at least one depression in the X-ray image within 
the inner regions (see Fig.[2])- Initially, we identified two cav- 
ities, one to the north-west and another to the south-east. 
However, as mentioned in Section 3, the BCG is located al- 
most within the south-eastern cavity. This could be due to 
projection effects, but we choose to discard this cavity and 
only consider the north-western one. If the cavity to the 
north-west is real, then the power capability is quite large 
(i^buoy ^ 4 X 10^^ erg s""*^ ). The optical image also shows 
an interesting structure with a central galaxy hosting 2 or 3 
potential nuclei. 

8.2.7 MACS J2245. 0+2637 

Finally, we present the results for MACS J2245.0+2637. The 
HST image shows a central dominant galaxy, with a very 
bright core. This cluster has a radio source associated with 
the central regions, as well as a short central cooling time of 
^ 4 Gyrs (Lcooi ^ 6 x 10^^ erg s~"^ ). We find one cavity in 
this system, but the data quality remains poor. The cavity is 
also quite small (^8 kpc in size) and the power stored within 
it is only Pbuoy 1 X lO^'^ erg s"^ or Pc^ 2 x 10^^ erg s"^ . 
This clu ster is also part of the 34 most luminous MACS 
clusters (jEbeling et al.ll2010h . 



9 DISCUSSION 
9.1 Selection effects 

We have analysed the Chandra X-ray observations of the 
MACS cluster sample and searched surface-brightness de- 
pressions in these systems. The MACS cluster sample con- 
sists of luminous, massive clusters all above a redshift of 
z > 0.3. Out of the 76 clusters, we find that 13 have "clear" 
surface-brightness depressions associated with their central 
AGN. We interpret these depressions as cavities being in- 
flated by the central SMBH, which then rise buoyantly into 
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Figure 10. Shown is the coohng luminosity of the system as a 
function of the average cavity radius in each system (if there are 
two cavities in one system, then we calculate and plot the average 
radius between the two cavities). The light blue po ints are for the 
MAC S clusters and the black points are for the iRaffertv et al.l 
(|2006i ) objects. 



the outskirts. If we include the 7 systems in which we find 
"potential" surface-brightness depressions, then the detec- 
tion rate goes up to 26 per cent. We stress that many of 
the clusters in which we identified cavities were those that 
had the "deepest" Chandra observations (more than 30 ks), 
and that many of the remaining MACS clusters only have 
10 — 20 ks observations (or as little as 500-1000 counts within 
a radius of 200 kpc). Our detection rates should therefore 
be regarded as lower limits. 

The 76 clusters that have been observed with Chandra 
were mostly selected on a flux limited basis te beling et al.l 
l2007l . FolO), and should therefore not be more biased to- 
wards cool cores than non cool cores. As shown in Fig. |4l we 
find the same percentage of cool core clusters in the MACS 
sample (with Chandra data) than non cool core clusters. 
However, the bottom panel of Fig. [5] shows that most of 
the clusters with deep observations are cool core clusters. 
This is because many of those that were initially observed 
for short exposure times and that showed interesting fea- 
tures, such as a peaked X-ray distribution, were followed- 
up with longer observations. We could therefore be missing 
many cavities in the non cool core systems which do not 
have long observations. However , as shown by earlier studies 
dBirzan et al."2004'; 'Pun n et al.1 12005': 'Raffertv et aD l2006l : 
iDunn & Fabian 2006, 200i), X-ray cavities seem to lie pre- 
dominantly in systems that require some form of heating, i.e. 
cool core clusters. Note also, the Chandra observations tar- 
geted the brightest clusters in the sample. We might there- 
fore be missing many of the less energetic outflows found in 
the less luminous clusters. This does not affect our results, 
since we draw our conclusions only on the powerful AGN 
outbursts and compare these to the powerful outbursts at 



lower redshift, both of which have Ln 



10 



'^^-^^ erg s-^ 



Finally, we mention that although the MACS sample 
lies within 0.3 ^ z ^ 0.7, the number of clusters with 
Chandra data decreases steadily with redshift and there are 
only 4 MACS clusters with Chandra data beyond z > 0.55. 
The probability of detecting cavities in systems beyond 



z = 0.55 is therefore low considering the limited sample. 
According to our cosmology, 1^^ corresponds to 4.45 kpc at 
z = 0.3 and 7.15 kpc at z = 0.7. Cavities have typical sizes 
of 20 kpc in radius (see Fig. [8]), whereas Chandra has a 
pixel size of ~ 0.49^^ and a point spread function of ^ l^^ A 
cavity of 20 kpc in radius located at z = 0.3 would there- 
fore spread across 4.5^^ (or 9 pixels), and if the same cavity 
were located at 2; = 0.7, it would spread across 3^^ (or 6 pix- 
els). Although deeper observations would be needed, and we 
would be pushing the unique resolution of Chandra to its 
limit, we could in principle identify cavities in ^ = 0. 6 — 0.7 
clusters, and might therefore be missing cavities in these 
further systems. 



9.2 Detection rates 

In more recent years, there has been considerable in- 
terest in studying the correlation between jet power, as 
determine d from cavity en e rgetic s , and radio luminosi- 
ties (e.g. [Birzan et al.1 12OO4I . I2OO8I : Icavagnolo et al.1 12OIOI : 
lO'Sullivan et al.ll201lL and see Fig.[7]here). This relation ex- 
tends from individual elliptical galaxies to groups and even 
clusters of galaxies, and is interesting for two reason. First, 
it pro vides insight into the nature of jets (e.g. IWillott et al.l 
"1999^), and secondly, it provides a means to determine jet 
powers from the rad^o properties alone fe.g. lBest et al.|[2007l : 
Magliocchetti & Briiggenl [20071) . Although the scatter re- 
mains large and more effort is needed to fully understand its 
nature, this method could in principle be used to determine 
jet powers in more distant clusters where the X-ray image 
quality remainsJx)o poor to identify surface-brightness de- 
pressions fsee lMa et al.ll201lL who looked at AGN heating 
in the 400SD survey using this relation). However, deriv- 
ing jet powers directly from X-ray cavities remains the most 
direct and reliabl e method. 

Birzan et all (l2004l) initially compiled a study of 80 
systems with Chandra X-ray archival data, in which the 
authors visually searched for surface-brightness depressions 
associated with the central regions. They found 18 sys- 
tems with well-defined cavities (16 clusters of galaxies; one 
galaxy group, HCG 62; and one giant elliptical, M84), cor- 
responding to a dete ction rate on the order of 20 per cent. 
iRaffertv etlH (|2006l ) then followed-up this study, by adding 
13 more objects where cavities had been recently discov- 
ered. Although many targets overlap between the various 
samples, Dunn Fabian (2006.) also provided an in depth 
study of X-ray cavities in t he Brightest 55 sample of galaxy 
clusters tedge et al.lll990l \ which was then followed by a 
more extensi ve study in clu ding th e Brightest Cluster Sam- 
ple (Ebeli ng et"al .11 19981) inlPunn &; FabianI (|2008i ) . The to- 
tal sample of IPunn k FabianI (|2006l . l2008h consists of 71 
clusters in the redshift range ^ z ^ 0.4. The authors 
find 22 systems with clear cavities, similar to our detection 
rate if we include the "potential" cavities (^ 25 per cent). 
However, Dunn FabianI (|2006l ) concluded that cavities sit 
mostly in clusters which require some form of heating, i.e. 
cool core clusters. They find that at least 14 of 20 (70 per 
cent) cool core clusters have cavities, making the detection 
rate in these systems very high. Since the majority of their 
clusters have an X-ray luminosity above 10^^ erg s~^ , simi- 
lar to the MACS clusters, we can directly compare our rates 
with theirs. We find that 20 (26 per cent), or at least 13 (17 
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per cent) of the 76 MACS clusters have cavities. We also find 
that the majority of our systems with cavities have a short 
central cooling time (see Fig. [3] and Fig. [4|. If we consider 
only the 37 MACS clusters that require some form of heating 
with tcooi(r<50kpc) < 3 Gyrs, then we find that 19 have cav- 
ities. The detection rate is therefore 50 per cent in cooling 
systems, but could be higher considering that more than half 
of the remaining 18 MACS clusters with tcooi(r<50kpc) < 3 
Gyrs have less than 30 ks observations (Fig. [5j). Note also, 
cool core clusters are more likely to be the ones with bright 
X-ray cores, since the central densities are high. When a cav- 
ity is present, the contrast it creates in these peaked cores is 
higher, which makes the identification of cavities easier. We 
might therefore be missing cavities in non cool core systems. 

More recently, studies have focussed on extending our 
knowledge of AGN feedback to the general population of 



MACS cluster sample 



giant ellipticals (e.g. Finoguenov & Jones 2001; Jone s et al 
2002; Machacek et al. 2006; Nulsen et al. 2007; Dunn et al 



I2OIOI ). iNulsen et al.l ()2007l ) searched through an extensive 
sample of 160 nearby giant ellipticals, and found 109 with 
significant emission from surrounding hot gas. From these, 
they find that 27 have X-ray cavities. This corresponds to a 
detection rate of only 25 per cent in these cooling systems 
(compared to 70 per cent in cool core clusters), suggesting 
that duty cycles are greater in larger systems (see Jones et 
al., in preparation). 



9.3 Cavity properties 

In the top panel of Fig. 2] we show that most of the sys- 
tems where we find cavities have a strong cool core with 
tcooi(r=50kpc) < 3 — 5 Gyrs, therefore requiring extreme feed- 
back from their central AGN to prevent the gas from cool- 
ing (Lmech ]^q44-45 s""*^ ). Although we expect to see 
some of the most extreme AGN outbursts in this sample, 
the cavities we find are not larger or more powerful than the 
counterparts found in nearby [z < 0.3) clusters with similar 

Fig. [6] shows that all but onj^ of the MACS clusters in 
which we find cavities are consistent with the Lmech — ^cooi 
correlation, implying that highly-luminous clusters are able 
to prevent the gas from cooling at the same l evel than the 
less- luminous clusters and groups. By contrast. iNulsen et al.l 
(I2OO7) provided evidence of a turnover at the lower lumi- 
nosity end, with outbursts being more powerful for a given 
cooling luminosity. This is consistent with the idea that in 
smaller systems, the duty cycles are smaller (only ^ 25 
per cent of the time, compared to ^ 70 per cent in cool 
core clusters), requiri ng that the outbursts be more ener- 
getic for a given cycle. iDunn k, FabianI (|2008h looked at the 
cavity properties of 6 distant clusters (0.1 ^ z ^ 0.4) and 
found that they all sit above the line of exact balance in the 
^mech — i^cooi corrclatlou. Our more complete study of dis- 
tant clusters finds differently. Essentially, our clusters agree 
well with the Lmech — ^cooi correlation with some above and 



2 The outlier refers to MACS J1359. 8+6231, where we found a 
cavity to the north-east, although very small (some 3 — 4 kpc in 
size) and therefore not energetically capable of preventing the gas 
from cooling. 
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Figure 11. In the top panel, we show the distribution of clusters 
in the MACS sample as a function of redshift. We outline those 
in which we identified cavities (also shown in more detail in the 
bottom panel). The fraction of clusters with cavities for each red- 
shift bin is shown in red. The fraction remains around 25 — 30 per 
cent until a redshift of 0.55. 



some below the line of exact balance. There is therefore no 
turnover at the higher luminosity end. 

In Fig. [HI we show the distribution of mechanical power 
(^mech), cavity energetics (pVtot) and buoyancy rise time 
(^buoy), as well as the cavity radius (i?average) and distance 
from the nucleus (i^distance)- According to these figures, we 
find that the MACS clu ster outbu rsts follow the sa me trends 
as the clusters in Raffertv et al.l (|200 6). The R affertv et alj 
(|2006l ) sample consists of clusters that are z ^ 0.3, some of 
which require the same heating powers as the MACS clusters 
to prevent the gas from cooling (Lmech ^ iq^^-^^ Q^g s~"^ ). 
In other words, the MACS outbursts are powerful, but not 
more powerful than the lower redshift counterparts. 

It is also important to emphasize that we can only in- 
terpret the results concerning the largest and most powerful 
outbursts. We are probably missing a number of smaller cav- 
ities which could not be detected since our clusters are at a 
relatively high redshift and the resolution of Chandra can 
only go so far. We will have to wait for future missions with 
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better spatial resolution than Chandra to study the small 
cavities at high redshift. 

Fig. [8] also shows that the pVtot distribution rises 
smoothly, peaks at 10^^~^° erg and then decreases steeply. 
By contrast, Lmech seems to rise and decrease more 
smoothly, although this distribution would be altered if 
we included the giant ellipticals for which Lmech peaks at 
-|^q4i-42 g-i Q^iYiQQ pVtot depends only on the thermal 
pressure and cavity size, the sudden drop at IQ^^"^'-' erg 
might be telling us something about how cavities form. As 
they are being inflated by the central AGN, it could be that 
they reach a maximum size in which they no longer can ex- 
pand no matter the physical conditions. There could there- 
fore be a maximum pVtot value that can be reached in a 
single outburst. To provide further insight into this idea, 
we plot in the left panel of Fig. [9] the energy in each cav- 
ity (pVtot) as a function of the average radius (-Raverage — 
(i?w X R\)^'^). Since pVtot oc i^average^, wc cxpcct a Correla- 
tion between these two quantities, and in a log- log plot, the 
slope should be 3 (this is shown with the dotted-line). How- 
ever, this plot also shows that the most powerful outflows 
can only be created by the largest cavities. If a small cavity 
was being created in a high-pressure environment, the result- 
ing pVtot could be large. However, Fig. [9] shows that we do 
not observe such cavities. Furthermore, the middle and right 
panels of Fig. [9] show that powerful outbursts are created in 
a variety of pressure environments, and that the size of each 
cavity has no dependence on the surrounding pressure. In 
other words, large cavities are created both in high and low 
pressure environments and the maximum allowed cavity ra- 
dius in clusters does not depend on the surrounding pressure. 
By contrast, one could argue that large cavities can only be 
created if the surrounding thermal pressure is not large. Fig. 
[9] suggests otherwise. Finally, in Fig. IIQI we plot the cooling 
luminos ity as a function of the average cavity radii in each 
system. iDunn FabianI (|2008) provided some evidence for 
a trend between these two quantities, and although the scat- 
ter remains large, our data support this idea (the larger the 
X-ray cooling, the larger the cavities). Again, this empha- 
sised that somehow the central AGN knows about the X-ray 
cooling and provides the cavities needed to offset cooling of 
the ICM down to a fine level. 

As shown by Fig. [H the distribution of cavity radii and 
distances from the nucleus also suggest that extremely large 
and distant cavities, suc h as those found in MS 0735.6+7421 
(|McN amar a et al.1 l2005l: iGitti et al. |i2007<l. are rare (these 
cavities are the outliers seen in the /^average and /^distance dis- 
tributions of Fig. [8]). Furthermore, /^distance /i^average peaks 
strongly at a value of 2 and then falls off quickly. In other 
words, cavities travel up to a distance of 2xi^average (i-e. 
their diameter) and then eithe r disintegjrate or, more likely , 
become too difficult to detect (jMcNamara Nulsenll2007h . 
More distant cavities have been detected in the Perseus 
cluster, b ut only through v ery deep X-ray observations 
(- 500 ks: ISanders FabianI E)07l : iFabian et al.ll201lh . 

Finally, we mention that the shapes of the cavities we 
find vary from almost perfectly circular to elongated and 
even flattened (e.g. MACS J2046.0-3430). Two clusters also 
seem to have a misalignment between the X-ray peak and 
central galaxy (MACS J2046.0-3430, MACS J2135. 2-0102), 
which may simply be due to projection effects. If a cavity is 
viewed along a different line of sight, projection effects can 



alter its shape and make it difficult to detect. The cavities we 
find may therefore be much larg:er and even mor e numerous 
than what we find (see also lSanders et al"]l2009l ). There are 
also two MACS clusters with long observations (> 100 ks) 
and short central cooling times (tcooi(^ < 50kpc) < 3 
Gyrs), but no cavities (MACS J131 1.0-0311 at z = 0.49 and 
MACS J1621.3+3810 at z = 0.465; see top panel in Fig. 
[5]). Statistically, we expect that some systems will have a 
jet axis aligned with the line of sight, rendering the cavities 
almost impossible to detect. This could explain why we do 
not see cavities in these two systems. 

9.4 Evolution of cavity properties 

The MACS sample consists of clusters within 0.3 ^ ^ ^ 0.7, 
for which we were able to find surface-brightness depres- 
sions in objects up to ^ ^ 0.6. In the top panel of Fig. [TT] 
we show the fraction of MACS clusters with cavities as a 
function of redshift, and in the bottom panel we show the 
redshift distribution of those with cavities. This fraction re- 
mains constant around 0.25 — 0.30 up to z ^ 0.6, and is 
consistent with the fra ction found in previous cavit y studies 
bv Birzan et al.l (|2004l . ^ 30 per cent) and R affertv et al.l 
(|2006l . ~ 25 per cent). Therefore, we find no significant evo- 
lution in the fraction of clusters with cavities as a function 
of redshift. 

An interesting finding in this study has been the lack 
of any significant evolution in the cavity properties. This is 
shown in Fig. 1121 where we plot the mechanical and energy 
outputs (Lmech and pVtot), as well as average cavity radius 
(i^average) and dlstaucc from the nucleus (/^distance) for each 
cavity, as a function of redshift. This figure shows that at 
least for the largest and most powerful outbursts, there is 
no evidence for evolution in the cavity properties. These 
cavities are not larger (or smaller) at higher redshift, and 
the distance to the nucleus remains on average the same, i.e. 
cavities at higher redshifts are not able to rise more (or less) 
than their lower redshift counterparts. The cavity energetics 
of the most powerful outbursts also remain the same. 

Therefore, AGNs lying at the centres of clusters have 
not only been able to power extreme mechanical outbursts 
as early as ^ 7 — 8 Gyrs after the Big Bang (correspond- 
ing to a redshift of z ^ 0.6), but they have also been 
able to maintain these extreme outbursts up to present day 
(Lmech(max) ^ 10^^ crg s~"^ ). lu a cosmological context, this 
means that "radio mode" feedback has been in place for 
at least the past 5 Gyrs and shows no sign of evolution 
since then. B ased on the CDM cosmological simulations of 
ICroton et al ] |2006), "radio mode" feedback occurs at late 
times and involves a phase where the AGN is accreting at 
sub-Eddington rates and driving powerful outflows that sup- 
press cooling flows from z ^ 1 to present. In this phase, it is 
thought that black holes do not grow significantly in mass. 
However, in order to maintain the powerful MACS outbursts 
for the past 5 Gyrs, the black hole mass must have been ini- 
tially very high. If we equate the energy required to power 
the outflows over the last 5 Gyrs (assuming an average out- 
put of ^ 10^^ erg s~^ ), to the energy released by a black hole 
(?7Mbhc^, assuming an efficiency of 77 = 0.1) then it takes 
at least a 10^ M© black hole to power such outflows. Here 
we assume that the duty cycles of dominant cluster galaxies 
are high (> 70 — 90 per cent), as supported by the observa- 



© 2011 RAS, MNRAS 000, ??-[2T] 



20 J. Hlavacek-Larrondo, et al. 



10000 



1000 



100 



10 



r 


1 1 1 1 

* ★ 


• • 












r •* 


* * 


• 








f 


* MACS: pot. bubbles 




MACS: cleor bubbles^! 




Rofferty et ol. 2006 
1 1 1 1 ■ 



0,0 0.1 0.2 0.3 0.4 0.5 0.6 

Redshift (z) 




0,0 0.1 0.2 0.3 0.4 0.5 0.6 

Redshift (z) 




0,0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 

Redshift (z) Redshift (z) 



Figure 12. Profiles of the mechanical energy being injected by the central AGN (top-left), the total energy output stored in the cavities 
(top-right), the average radius of a bubble (bottom-left) and distance from the core to the cavity centre (bottom-right) as a function of 
redshift. The MACS cl u sters are shown with the star symbols (dark blue for "clear" cavities and light blue for "potential" cavities), and 
those of iRaffertv et al are shown with the grey circles. Objects common to both samples are shown in red and yellow. 



tion that most or even all of the clusters that require some 
form of heating to counterb alance cooling have jetted out- 
flows (jPunn k, Fabianll2006l ). Note also, at least 50 per cent 
of our clusters requiring feedback from their central AGN 
seem to have jetted outflows, therefore supporting the idea 
that duty cycles are high (see Section 9.2). Even with a duty 
cycle of 50 per cent, the black hole masses required to power 
the outflows are on the order of 10^ for an efficiency of 
Tf 0.1. This rivals the largest known black holes masses, 
and suggests that some, especially those lying at the centres 
of the most extreme cooling flows, may be very massive and 
even ultramassive with > lO^^M©. 

Likewise, we could argue that to power outflows of 
10^^ erg s~^ for the past 5 Gyrs, this requires the accreted 
mass (Mace = ^cav/'/7c^) to be at least on the order of 
lO^M©. If the initial mass of the black hole was ^ 10^ Mo, 
then by accreting some 10^ M© in mass, we should have seen 
some variation in the cavity properties with redshift. Indeed, 
the mass of the black hole is one of the fundamental prop- 
erties that regulates the power capabilities (another would 
be the spin). Since we do not see any evidence for evolu- 
tion in the outflow properties, this strongly suggests that 
the 10^ Mo increase in mass does not signiflcantly affect the 



AGN, and therefore that the black holes are initially very 
mas sive with Mbh > lO^^^M©. 

ISamuele et all (|201lh looked at nebular emission in 
77 BCGs belonging to the 160SD X-ray cluster survey 
(IVikhlinin et al.M 19981 : iMullis et al] 2003) and found no 
[Oii]A3727 or Ha emission stronger than 15A and 5 A re- 
spectively in any of the BCGs. The authors suggest that this 
may be due to a signiflcant decrease in the number of strong 
cooling flow clusters from z ^ 0.5 to today (which does not 
mean that there are no strong cool cores, just less). They 
also suggest that this decline could be due to over-hea ting 
of the core by AGN (as supported by i Ma et al.1 l201lh or 
cluster mergers, which happen more frequently at higher 
redshift. By using the relation between jet power and ra- 
dio power (Birzan et al. 2004, 2008; Cavagnolo et al. 2 oTol : 
O'Sul livan et al.ll201lhJMa et al.1 (|201lh find that the power 
injected by radio sources within r < 250 kpc is significant 
compared to the power radiated by the cluster atmosphere. 
However, as shown in Fig. [6l at least for the central AGN 
in the MACS clusters, the outbursts we find are able to bal- 
ance cooling down to a fine level up to z ^ 0.6. The central 
AGN in clusters therefore does not seem to overheat the core 
at higher redshift and can not account for the decrease in 
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strong cool cores with redshift. Other studies in the htera- 
ture also find a decrease i n the number of strong cooling flow 
clusters with redshift fsee lVikhlinin et al ] l2007l :l ISantos et alJ 
'2008), but suggest that it is most likely due to the higher 
major merger rate of the past. 

As a last note, we mention that many of the clusters in 
which we find cavities have a bright X-ray nucleus associated 
with the central AGN. These include MACS J0547.0- 
3904, MACS J0913.7+4056, MACS J0947.2+7623, 
MACS J1411.3+5212, MACS J1423.8+2404, 
MACS J1931.8-2634 and MACS J2046.0-3430. Most 
of these either have an Fe K line in the spectrum of 
the nucleus or a spectrum that can best be fitted by a 
power- law, indicating the presence of an active AGN. BCGs 
with bright X-ray cores are rare at low redshift, and many 
that have similar Lmech as those in the MACS clusters 
do not have any evidence of an X-ray nucleus associated 
with their central AGN (Hlavacek-Larrondo & Fabia^ 
l201lh . This may be indicating that we are seeing some 
evolution, not in terms of the mechanical properties of AGN 
feedback, but in terms of the radiative properties of the 
central AGN, thus suggesting that we are starting to see 
the transition between "quasar mode" and "radio mode" 
feedback. This will be analysed in detail in a forthcoming 
paper (Hlavacek-Larrondo et al., in preparation). 



10 CONCLUDING REMARKS 

We have analysed the Chandra X-ray observations of the 
MACS cluster sample and searched for surface-brightness 
depressions associated with the central AGN in each cluster. 
The MACS sample consists of very X-ray luminous clusters 
within 0.3 ^ z ^ 0.7 (median Lx,rass = 7 x 10^^ erg s~^ ). 
Out of the 76 clusters, we find 13 with "clear" cavities and 7 
with "potential" cavities, bringing the detection rate to ^ 25 
per cent. Most of the clusters in which we find cavities have a 
short central cooling time below 3 — 5 Gyrs, consistent with 
the idea that cavities sit predominantly in cool core clus- 
ters. By combining our results with those of previous cavity 
studies, the latter focussing on systems at z < 0.3, we find 
no evidence for evolution in any of the cavity properties. 
Although we expect to see some of the most extreme out- 
bursts (Lmech ^ io^^~^^ erg s~"^ ), the cavities are not larger 
(or smaller) at higher redshift, and they are not able to rise 
to further (or lesser) radii. The energy and power capabil- 
ities of the most powerful outbursts also remain the same. 
Extreme "radio mode" AGN feedback therefore starts to op- 
erate as early as 7 — 8 Gyrs after the Big Bang and shows 
no sign of evolution since then. In other words, AGNs lying 
at the centre of clusters are able to operate at early times 
and with extreme powers. 
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APPENDIX A: UNSHARP-MASKED AND 
ELLIPSE-SUBTRACTED X-RAY IMAGES 

We show the unsharp-masked and ellipse-subtracted images 
for all clusters in which we identified cavities fFig. lAl)) . The 
first consists of subtracting a strongly smoothed image from 
a lightly smoothed image, and the second consists of sub- 
tracting an elliptical model of the cluster emission from the 
original image. These techniques enhance deviations in the 
original image, and the images were used as a first step to 
identify cavities. However, we only considered that a cluster 
had a cavity if we could also see a faint depression in the 
original 0.5 — 7 keV Chandra image. There are two clus- 
ters for each horizontal line in Fig. IA1[ with the unsharp- 
masked 0.5 — 7 keV X-ray image shown in the left panel and 
the ellipse-subtracted 0.5 — 7 keV X-ray image shown in the 
right panel. The red cross shows the location of the BCG 
nucleus. We also indicate in the lower-left corner of each 
image the binning (Bx) and smoothing (Sx) scales adopted 
for creating the images. For Bx, X indicates the binning 
factor where X = 1 corresponds to no binning and X = 2 
gives an image where each of pixel corresponds to 4 pixels 
in the original image. For Sx, X indicates the smoothing 
factor and corresponds to the sigma of a gaussian in units 
of pixels. For each unsharp-masked image, we show the 2 
smoothing scales used to create the image (Sxi — Sxs), and 
for each ellipse-subtracted image, we show the smoothing 
scale adopted (Sxs) for the original image before creating 
and subtracting an elliptical model. The first 13 clusters 
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Figure Al. See caption on next page. 
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Figure Al. Unsharp-masked and ellipse-subtracted images for all clusters in which we identified cavities. Same layout and scaling as Fig. 
[21 two clusters for each horizontal line with the unsharp-masked 0.5 — 7 keV X-ray image shown in the left panel and the ellipse-subtracted 
0.5 — 7 keV X-ray image shown in the right panel. We identify the cavities in each system, and we show with the red cross symbol in the 
X-ray image, the location of the BCG. We also show in the lower-left corner of each image the binning (Bx) and smoothing (Sx) scales 
adopted for creating the images. For each unsharp-masked image, we show the 2 smoothing scales used to create the image (Sxi — 8x2)7 
and for each ellipse-subtracted image, we show the smoothing scale adopted (8x3) for the original image before creating and subtracting 
an elliptical model. If no Sx is shown, then no smoothing was applied to the image. The first 13 clusters have "clear" cavities, while the 
remaining 7 only have "potential" cavities (see "potential cavities" indicated in the figure). We emphasize that deeper observations are 
needed to confirm if the "potential" cavities are real. 



have "clear" cavities, while the remaining 7 only have "po- 
tential" cavities. Deeper observations are needed to confirm 
if the "potential" cavities are real. 



APPENDIX B: SPECTRAL PROFILES 

Fig. lBll shows the projected (black symbols) and deprojected 
(red symbols) temperature, electron density, electron pres- 
sure and entropy profiles for each cluster in which a cavity 
was identified. We used statistics and each annulus had 



a minimum of 900 counts, but for clusters in which the data 
was of good quality we allowed the annuli to have up to 4000 
counts. For MACS JOlll. 5+0855 and MACS J2135. 2-0102, 
we could not deproject the spectra because of the limiting 
data quality. In this case, we used C-statistics and each an- 
nulus only had ^ 400 counts. 
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Figure Bl. Continued - Projected (black symbols) and deprojected (red symbols) temperature, electron density, electron pressure and 
entropy profiles for each cluster in which a cavity was identified. 
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